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PREFACE 


The  ultrasonic  harmonic  generation  technique  has  been  extended 
to  the  study  of  the  nonlinear  elastic  behavior  and  thermal  properties 
of  the  diamond-like  solids  silicon  and  germanium.  A  complete  study  of 
the  temperature  variation  of  the  third-order  elastic  constants  of 
silicon  and  germanium  has  been  done.  In  order  to  make  the  technical 
report  intelligible,  it  has  been  divided  into  four  chapters  as  follows. 

The  importance  of  the  study  of  the  nonlinear  properties  of 
solids  is  emphasized  in  the  Introduction  which  forms  the  first  section 
of  Chapter  I.  A  brief  review  of  the  theoretical  and  experimental  work 
done  in  higher-order  elasticity  is  given.  The  remainder  of  Chapter  1 
concentrates  on  the  theoretical  basis  of  the  work.  The  nonlinear  wave 
equation  in  a  cubic  solid  is  derived  and  its  solutions  alonq  the 
symmetry  directions  are  given.  The  ultrasonic  nonlinearity  parameter 
for  lonqitudinal  wave  propagation  along  pure  mode  directions  in  cubic 
crystals  is  defined.  The  relationship  between  the  nonlinearity 
parameters  and  elastic  constants  is  derived. 

Chapter  II  describes  the  experimental  setup.  A  brief  description 
of  the  capacitive  detector  which  is  used  for  wave  amplitude  measurements 
is  given.  Both  room  temperature  and  cryogenic  apparatus  are  described 
with  diagrams  wherever  necessary.  The  measurement  procedure  is  also 
outlined  in  this  chapter.  The  block  diagrams  of  the  room  temperature 
as  well  as  cryogenic  apparatus  are  given.  Details  about  the  individual 
parts  are  qiven  wherever  necessary.  Ultrasonic  nonlinearity  narameters 
of  silicon  have  been  measured  at  room  temperature  and  as  a  function  of 
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temperature  down  to  liquid  helium  temperature.  The  numerical  results 
are  tabulated  and  the  temperature  variation  of  the  nonlinearity 
parameters  and  some  combinations  of  third-order  elastic  constants  have 
been  plotted.  Chapter  III  is  made  up  of  these  results.  Comparison  of 
the  results  with  previous  experimental  and  theoretical  values  are  made 
as  far  as  possible.  The  strain  generalized  Griineisen  parameters  along 
the  principal  directions  are  evaluated  and  reported. 

Keating's  central  force  model  is  one  of  the  most  successful 
and  established  models  for  diamond-like  solids.  We  have  developed  a 
method  to  isolate  all  the  six  individual  third-order  elastic  constants 
from  nonlinearity  parameters  of  diamond-like  solids  bv  combining  our 
measured  values  with  this  lattice  dynamical  model.  Expressions  for  the 
Keating  force  constants  are  given  in  terms  of  the  nonlinearity  parameters. 
The  method  is  apolied  to  silicon  and  germanium  and  the  temperature  variation 
of  all  the  six  independent  third-order  elastic  constants  have  been 
studied  between  room  temperature  and  3°K.  The  method  and  results  are 
given  in  Chapter  IV.  Numerical  knowledge  of  the  temperature  dependent 
third-order  elastic  constants  can  be  related  to  the  temperatu.v  dependence 
of  the  anharmonic  parameters  such  as  thermal  expansion,  Griineisen 
parameter,  etc.  of  solids.  Such  a  study  made  for  silicon  and  germanium 
also  is  incorporated  in  this  chapter.  Low  temperature  thermal  expansion 
measurements  of  diamond-like  solids  by  various  workers  have  shown  that 
they  possess  a  negative  thermal  expansion  at  low  temperatures,  Gruneisen 
parameters  follow  the  same  anomaly.  Theoretical  calculations  of  low 
temperature  thermal  expansion  and  Gruneisen  parameters  from  room 
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temperature  third-order  elastic  constants  by  earlier  workers  do  not 
account  for  the  negative  thermal  expansion  in  certain  temperature  ranqes. 
We  have  tried  to  bridge  this  gap  by  evaluatina  the  thermal  expansion  and 
Griineisen  parameters  of  silicon  and  germanium  as  a  function  of  temperature 
from  temperature  dependent  third-order  elastic  constants  in  the  ouasi- 
harmonic  approximation.  Temperature  variation  of  other  anharmonic 
parameters  of  silicon  and  germanium  have  been  studied  and  the  results 
are  presented. 


CHAPTER  I 


INTRODUCTION  AND  THEORY 

A.  INTRODUCTION 

A  solid  medium,  in  general,  is  nonlinear.  Hooke's  law  is  not 
obeyed  by  solids  because  of  this  nonlinearity.  Nonlinearitv  arises 
from  the  hiqher  order  terms  in  the  eneray  density  expansion  of  the 
solid.  A  finite  amplitude  ultrasonic  wave  nropaqating  throunh  a  non¬ 
linear  medium  distorts  as  it  progresses  and  Generates  hiqher  harmonics. 
Superposed  waves  qenerate  sum  and  difference  frequency  waves.  The 
nonlinear  effects  are  mainly  due  to  the  asymmetry  of  the  interatomic 
potential  well;  i.e.,  due  to  a  nonparabol ic  interatomic  interaction 
potential  which  leads  to  anharmonic  behavior  of  the  crystal  lattice. 

The  distortion  of  a  finite  amplitude  ultrasonic  wave  in  a  solid  is  a 
measure  of  the  nonlinearity  of  the  medium  and  measurement  of  the 
Generated  harmonics  offers  a  unique  method  for  the  study  of  anharmonici tv 
of  solids.  Distortion  of  a  finite-amplitude  ultrasonic  wave  in  solids 
was  first  observed  by  Breazeale  and  Thompson^  and  by  Gedroits  and 

p 

Krasilnikov.  '  The  fact  that  such  measurements  lead  to  the  determination 
of  third-order  elastic  constants  has  given  impetus  for  the  work  to  be 
described  in  the  succeeding  chapters. 

Tn  describing  the  physical  properties  of  solids,  the  elastic 
constants  play  a  central  role.  In  the  macroscopic  theory  of  elasticity 
the  solid  is  treated  as  a  continuum  which  has,  in  the  unstrained 


equilibrium  state,  a  constant  mass  density.  I  he  elastic  enemy  density 

is  written  as  a  power  series  in  the  elements  of  the  strain  tensor,  lhe 

zero-order  term  represents  the  enerqy  of  the  solid  in  eouil ibrium.  It 

is  independent  of  the  strain  and  may  be  set  equal  to  zero.  The  first- 

order  term  is  also  zero  in  the  absence  of  external  body  forces.  The 

coefficients  of  the  second-,  third-,  and  hioher-order  terms  are 

designated  as  elastic  constants  of  order  correspondina  to  the  power  of 

the  strain  in  the  respective  terms.  In  the  microscopic  theory  of 

crystal  lattices,  the  potential  enerqy  of  a  crystal  is  expanded  as  a 

Taylor  series  in  terms  of  the  displacements  of  the  atoms  in  the  lattice 
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from  their  equilibrium  positions.  The  coefficients  in  the  series  are 
the  derivatives  of  the  potential  enerqy  referred  to  the  equilibrium 
positions  and  are  called  the  coupling  parameters.  The  zero-order  term 
may  be  set  equal  to  zero;  the  first-order  term  vanishes  because  the 
expansion  is  about  the  equilibrium  positions.  As  in  the  macroscopic 
theory,  the  higher  terms  stand  for  the  corresnondino  couplina  parameters. 
The  symmetrv  properties  of  the  lattice  can  be  used  to  obtain  the 
relations  amonq  the  coupling  parameters  and  the  corresponding  elastic 
constants.  Measurement  of  the  elastic  constants  yields  information 
about  the  lattice  forces. 

If  only  the  second-order  terms  are  retained  in  the  elastic  energy 
density  and  linear  elasticity  theory  is  used,  a  linear  relationship 
between  stress  and  strain  results  (Hooke's  law).  This  approximation 
can  be  used  for  infinitesimal  deformations;  i.e.,  when  the  snace 
derivatives  of  the  displacement  vectors  of  a  point  in  the  bodv  are  small 
enounh  that  their  products  and  squares  may  be  nealected  in  comparison 


to  the  terms  themselves.  Similarly,  neglectinn  terms  beyond  second 

order  in  the  lattice  potential  enerqy  leads  to  lattice  dynamics  in 
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the  harmonic  approximation  with  which  one  can  account  for  much  of  the 
behavior  of  solids. 

A  number  of  properties  of  solids  are  related  to  the  anharmonici tv 
of  the  lattice  potential  which  involves  elastic  constants  hioher  than 
second  order.  Crystals  exhibit  lattice  thermal  expansion.  The 
adiabatic  and  isothermal  elastic  constants  are ,  in  oeneral,  different 
from  each  other  and  the  elastic  constants  do  vary  with  temperature  and 
pressure.  At  high  temperatures  the  specific  heat  exhibits  a  temperature 
dependence.  All  these  anharnonic  properties  indicate  that  harmonic 
theory  is  not  sufficient  for  a  quantitative  description  of  the 
properties  of  a  crystal. 

Even  at  low  temperatures  the  harmonic  theory  is  not  correct. 
Because  of  zero-point  vibrations  the  influence  of  anharmonicity  does 
appear  even  at  absolute  zero.  This  is  particularly  large  for  small  atoms 
or  molecules  which  have  correspondingly  large  zero-point  oscillations. 

The  eouilibrium  positions  therefore  do  not  coincide  with  the  minimum 
of  the  potential  energy.  The  methods  of  harmonic  theory  are  no  longer 
applicable  in  the  dynamic  theory  of  elasticity  because  the  anharmonic 
terms  cause  a  coupling  between  the  different  modes  of  oscillations. 

In  the  harmonic  theory  the  oscillator  or  eigenfreauencies  wj  of  the 
lattice  are  independent  of  interatomic  distance  a.  With  the  correct 
potential  energy,  the  coupling  parameters  of  second  order  and  therefore 
are  functions  of  a.  Even  if  the  terms  higher  than  the  ouadratic 
ones  in  the  energy  expansion  are  omitted,  part  of  the  anharmonic  effects 
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are  already  described  by  the  dependence  oo,  (a).  This  approximation  is 

J 

called  the  quasi  harmonic  approximation  which  is  conveniently  used  for  the 
evaluation  of  many  physical  parameters.  Excellent  reviews  on  the  anhar- 
monic  effects  in  solids  are  available  in  literature.6 

The  conduction  of  heat  in  solids  depends  on  the  transfer  of 
vibrational  enerqy  through  the  lattice,  which  is  not  predicted  by  har¬ 
monic  theory.  Anharmonic  terms  represent  the  interaction  between 
mechanical  sound  waves  and  thermal  vibrations,  which  leads  to  the 

g 

damping  of  high  freauency  (>  10'  Hz)  sound  waves  and  to  the  damping  of 

moving  dislocations  by  phonon  viscosity.6  Volume  increase  of  an  elastic 

body  associated  with  a  change  of  temperature  or  the  introduction  of 

dislocations,  an  exact  description  of  stress  and  strain  fields  near 

dislocations,  the  electrical  resistance  of  screw  and  edge  dislocations, 

scattering  of  phonons  from  defects,  etc.  demand  a  treatment  by  nonlinear 
7  8 

elasticity  theory.  ’ 

The  elastic  constants  play  a  dual  role  of  importance  in  the 
theory  of  solids.  Only  a  knowledge  of  the  elastic  constants  enables  one 
to  determine  such  important  physical  properties  as  compressibility,  the 
Debye  temperature  and  its  pressure  dependence,  the  Griineisen  parameters 
and  various  thermodynamic  properties  of  defects  in  solids.  These 
properties  can  be  found  without  havinn  any  information  about  the  inter¬ 
atomic  forces  or  potentials.  The  second  important  role  of  elastic 
constants  is  a  result  of  their  being  among  the  most  accurately  measurable 
physical  quantities.  Ultrasonic  techniaues  have  generally  enabled  one 
to  experimentally  determine  elastic  constants  to  a  greater  degree  of 
accuracy  than  can  be  obtained  from  theoretical  calculations.  Thus,  the 
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elastic  constants  can  serve  as  a  useful  guide  in  de\  .oping  any 

related  theory.  Since  elastic  constants  describe  how  the  enerqy 

density  of  a  solid  changes  with  respect  to  various  volume  and  shear 

deformations,  they  are  useful  in  determinina  the  nature  of  bindino 

forces  and  interatomic  potentials.  The  macroscopic  nonlinear 

behavior  of  solids  is  determined  by  the  intermolecular  potential 

functions  of  solids.  The  relations  between  the  force  constants  and 

elastic  constants  of  a  cubic  solid  has  been  given  by  Coldwell- 
9  10 

Horsfall.  Ghate  uses  the  Born  model  of  ionic  solids  to  compute 
the  third-order  elastic  constants  and  their  temperature  variation. 

Hi ki  and  Granato^  investigate  short  range  repulsive  interactions  in 
noble  metals  and  their  relationships  to  the  nonlinearity  of  the  solid. 
Basically  there  are  two  methods  for  calculating  the  elastic 

constants  of  a  crystalline  medium:  the  method  of  homogeneous 

12  13 

deformation  and  the  method  of  lonq  waves.  The  results  obtained 

from  the  two  methods  should  be  in  agreement  when  the  same  model  of  a 

solid  is  used  in  both  cases.  In  the  method  of  homogeneous  deformation 

the  calculation  proceeds  by  determining  how  the  energy  density  changes 

with  respect  to  various  homogeneous  deformations;  i.e.,  deformations 

for  which  the  resulting  structure  remains  a  perfect  lattice.  The  two 

most  widelv  used  types  of  deformations  are  those  described  by  the 

Lagranqian  strain  parameters  and  by  the  Fuchs  strain  parameters.  In 

the  first  case,  derivatives  of  the  energy  density  are  taken  with 

respect  to  the  Lagrangian  strains  n...  The  resulting  derivatives  when 

*  «J 

evaluated  in  the  undeformed  state,  are  the  Bruqger  elastic  constants 


9 


Cijkt. . 


=  (:»nE/anij3nkjl...)n=0 


(1.1) 


where  E  is  the  energy  of  the  crystal  per  unit  undeformed  volume.  This 

thermodynamic  definition  of  elastic  constants  of  any  order  introduced 
14 

by  Brugger  is  most  widely  used  nowadavs.  In  the  other  case,  one 

calculates  the  derivatives  of  the  energy  density  with  respect  to  Fuch- 

1 5 

tvpe  strain  parameters  which  consist  of  a  volume  deformation  and  of 

various  volume  conserving  shear  deformations.  When  evaluated  in  the 

undeformed  state,  these  derivatives  are  referred  to  as  Fuchs  elastic 

constants.  The  Fuchs  constants  are  linear  combinations  of  the  Brugger 

constants.  The  relationships  between  the  Fuchs  and  Brugger  elastic 

16 

constants  for  cubic  crystals  are 

A/^2  =  C44 


A/iv2  =  1 

(cn  +  2Ci2) 

A/Sr..,2  =  2(CU  -  C 

12^ 

i3E/3V3  =  ^ 

(cm  + 

2C1 23 

+  6Cn2 

>  -  J 

(Cn  +  2C 

12^ 

1  (1 

^E/SVar-,2  = 

1  (cm 

-  r  )  +  — 
C123;  +  3 

(C11  " 

C12>  +  1 

(Cn  +  2c12) 

^3e/^v^12  = 

1  (C144 

+  2C 

166^  +  I 

(Cll 

+  2^1 2  ^  + 

4  r 

1  1 44 

33E/3f23Yi 2 

=  C1 66  ' 

C 1 44 

+  cn  - 

C1 2  + 

2C144 

33E/3£23f-  = 

1 

Clll  + 

2C1 23 

-  3C-|  1 2 

+  7  ( C 

11  "  C12^ 

1  *  C456 

. 

10 


'.j,  Y.j  and  V  are  the  Fuchs  strain  parameters.  Tor  most  symmetric  cubic 
crystals  there  are  three  independent  second-order  elastic  constants  and 
six  third-order  elastic  constants.  The  Voigt  notation  for  elastic 
constants  is  used  throughout  the  text  of  this  report.  The  relationship 
between  Fuchs  and  Rrugger  elastic  constants  for  hexagonal  crystals  have 
been  given  elsewhere.^ 

The  Fuchs  constants  are  generally  more  convenient  to  calculate 

when  dealing  with  noncentral  potentials,  whereas  the  Brugger  constants 

are  the  preferred  type  for  central  potentials.  In  either  case,  the  final 

results  are  usually  expressed  in  terms  of  the  Brugger  constants  to 

facilitate  comparison  with  experimental  data. 

The  elastic  strain  energy  for  most  symmetrical  class  of  cubic 

crystals,  including  third-order  terms,  but  omitting  terms  independent 

18a 

of  strain  was  first  given  by  Birch  as 


1  2  2  2 

<!>  =  o  cn(ri11  +  n22  +  n33 )  +  ci 2^nll n22  +  n22n33  +  n33nl  1  ^ 

.  r  ,  2  .  2  .  2,  2,  2  ,  2Wr  ,  3.  3.  3, 

C44'r,12  +  n21  n23  n32  n31  n13^  +  Clll'n11  n22  n33' 

+  ci  1 2^nl  1  ^n22  +  n33^  +  n22 ^ n33  +  nll  ^  +  n33^nll  +  n22^ 

+  1  C144^nl 1 ^n23  +  n32^  +  n22  ^n3?  +  n13^  + n33^n12  +  n21 ^ 

+  \  C166^n12  +  n21  ^nll  +  n22^  +  ^n23  +  n32^(n22  +  n33^ 

2  2 

+  (r'31  +  n13)(n33  +  )]  +  C123  n11n22n33 


+  C,icr[ 


456un12n23n31  +  n21n32n13^ 


(1.3) 


where  are  the  Lagrangian  strain  components  used  by  Birch  and  we 
have  corrected  the  numerical  errors  pointed  out  by  Bhagavantam  and 
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Suryanarayana .  The  relationships  between  Brugqer  and  Birch  third- 
order  elastic  (TOE)  constants  can  easily  be  worked  out  and  are  qiven  by 

C111  (Birch)  =  ^  C]11  (Brugger) 

(Birch)  =  ^  C-j^  (Brugaer) 

C144  (Birch)  =  2ci44  (Pruqger) 

(1.4) 

C166  (Birch)  =  2C16g  (Brugqer) 

C ^ ^ 3  (Birch)  =  (Brugger) 

C456  (Birch)  =  4C456  (Brugger)  . 

The  role  of  second-order  elastic  (SOE)  constants  in  the  description 

of  the  elastic  properties  of  solids  is  well  established  and  appears  in 

every  work  in  the  field  and  so  is  not  repeated  here.  The  nulse  echo 

technique  of  determining  the  velocity  of  ultrasonic  waves  in  crystals, 

1 9 

first  employed  by  McSkimin,  is  one  of  the  most  accurate  and  widely 
prevalent  methods  of  measuring  the  SOE  constants  of  a  cystal .  It 
involves  the  application  of  a  series  of  RF  pulses  of  about  15  to  25  MHz 
to  a  transducer  bonded  to  the  specimen  by  a  suitable  bond.  The  pulse 
repetition  rate  is  adjusted  so  that  the  pulse  echoes  superimpose  exactly 
on  the  siqnal.  This  pulse  repetition  frequency  (~  100  to  500  kHz)  is 
directly  proportional  to  the  square  of  the  natural  wave  velocity  and 

5 

can  be  measured  with  an  accuracy  of  1  in  10  .  Measurement  of  the 
velocities  of  longitudinal  and  shear  waves  along  the  symmetry  directions 
leads  directly  to  the  SOE  constants.  The  nulse  interference  method^  is 
also  often  used  for  velocity  measurements. 

In  explaining  the  anharmonic  properties  of  solids,  the  TOE 
constants  play  an  important  role.  They  are  necessary  to  evaluate  the 
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third-order  terms  of  the  lattice  potential  enerov  and  are  needed  to 
evaluate  the  generalized  Gruneisen  parameters  which  describe  the  strain 
dependence  of  the  lattice  vibrational  frequencies.  The  TOE  constants 
can  provide  quantitative  descriptions  of  acoustic  amplification  at  micro- 
wave  frequencies  besides  enabling  the  evaluation  of  the  Akhieser  phonon- 
phonon  interaction  mechanism  which  describes  the  attenuation  of  ultrasonic 
waves  in  solids.  The  TOE  constant  data  are  needed  to  determine  the 
changes  in  the  lattice  parameters  of  the  solid  due  to  application  of 
hydrastatic  pressure. 

The  most  common  method  of  determini  no  the  TOE  constants  is  to 

measure  the  changes  in  velocities  of  ultrasonic  waves  propaaatinq  alono 

20 

different  symmetry  directions  with  applied  stress.  Hughes  and  Kelley 

derived  expressions  for  the  velocities  of  elastic  waves  in  stressed 

solids  using  third-order  terms  in  the  elastic  enerqy  and  finite  strain 
18  21 

theory  of  Birch.  EinsDruah  and  Manning  applied  finite  strain 

elasticity  theory  to  evaluate  the  TOE  constants  of  anisotropic  crystalline 

solids  and  has  presented  results  for  cubic  and  uniaxial  crystals.  The 

first  measurement  of  TOE  constants  of  isotropic  materials  polystyrene, 

20 

iron  and  pyrex  glass  has  been  done  by  Hughes  and  Kelly  and  those  of  the 

22 

anisotropic  crystal  aermanium  by  Bateman  et  al.  The  evaluation  of  all 

the  independent  TOE  constants  involves  measurement  of  wave  velocities 

when  uniaxial  static  stress  is  applied  to  the  crystal,  besides  measure- 

23 

ments  under  hydrostatic  pressure.  Thurston  and  Brugoer  ’  have  derived 
the  expressions  for  the  sound  velocity  and  for  a  natural  velocity  and 
their  stress  derivatives  evaluated  at  zero  stress  in  terms  of  SOE  and  TOF 
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constants  of  the  crystal  applicable  to  arbitrary  crystal  symmetry  and 

arbitrary  stress  systems  dependina  on  a  single  scalar  variable. 

1 9 

In  the  pulse  echo  method  for  measuring  the  chanoes  in  wave 
velocity  due  to  application  of  hydrostatic  pressure,  the  specimen  is 
subjected  to  high  pressure  in  a  cylindrical  bomb,  using  comDre ssed 
gases  or  liquids.  Measurements  under  uniaxial  stress  are  performed  by 

3 

applying  the  stress  through  a  20  x  10  kg  hydraulic  jack  operated  by  a 

screw.  The  ratio  of  the  pulse  repetition  freauency  when  the  crystal 

is  stressed  to  the  repetition  frequency  of  the  unstressed  crystal  is 

plotted  against  the  stress.  From  the  slope  of  the  linear  portion  of 

the  graph  for  different  propagation  directions,  wave  polarizations  and 

stress  systems,  an  adequate  number  of  combinations  of  TOE  constants  can 

be  obtained.  A  set  of  simultaneous  equations  are  solved  to  obtain  the 

individual  TOE  constants.  The  expressions  connecting  TOE  constants  and 

stress  derivatives  of  sound  velocities  in  solids  have  been  aiven  by  a 
22  24  25 

number  of  authors  ’  ’  and  are  not  reproduced  here.  The  pressure 

technique  has  been  employed  by  a  number  of  later  workers  to  determine 

26  27  28 

the  TOE  constants  of  silicon,  alkali  halides,  fused  silica,  noble 

1 1  29 

metals,  etc.  Salama  and  Alers  used  the  technique  to  measure  TOE 

constants  of  copper  at  low  temperatures . 

30 

Graham  has  described  a  method  of  obtaining  longitudinal  TOE 
and  FOE  (fourth-order  elastic)  constants  of  solids  which  sustain  elastic 
compressions  under  shock  wave  loading.  The  shock  compression  data  on 
sapphire  and  fused  quartz  have  been  analyzed  to  determine  C^i,  C^, 
Gjlll  and  £3333-  The  determination  of  the  TOE  constants  under  these 
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larqe  compressions  allows  one  to  test  the  applicability  of  the  finite 

strain  formulation  of  constitutive  relations.  This  technique  can  be 

applied  to  materials  whose  elastic  limits  are  a  few  percent  of  their 

lonqitudinal  elastic  constants. 

An  optical  technique  for  determining  the  TOE  constants  of 

31 

transparent  crystals  has  been  proposed  and  measurements  have  been 
32 

carried  out  on  NaCl  and  C-j^  has  been  reported.  This  techniaue 

makes  use  of  the  fact  that  when  an  initially  sinusoidal  ultrasonic  wave 

of  finite  amplitude  is  propagating  in  a  crystal,  it  qets  distorted  due 

to  the  generation  of  harmonics.  The  difference  in  intensity  between 

the  first  positive  and  first  neqative  orders  of  the  diffraction  pattern 

obtained  when  monochromatic  light  is  transmitted  through  the  crystal 

perpendicular  to  he  direction  of  propagation  of  the  ultrasonic  wave  is 

a  function  of  the  strain  amplitude  of  the  wave.  A  measurement  of  the 

asymmetry  in  the  intensity  leads  to  some  of  the  TDE  constants.  Not  much 

work  has  been  done  in  this  regard. 

Measurements  of  phonon-phonon  scattering  have  been  used  to 

determine  various  combinations  of  TOE  constants.  Classical  calculations 

based  on  the  interaction  of  two  intersecting  ultrasonic  beams  within  an 

33 

isotropic  solid  have  been  carried  out  by  Jones  and  Kobett.  Ouantum 

mechanical  calculations  of  the  same  experimental  situation  have  been 

34 

treated  by  Taylor  and  Rollins  and  experimental  results  for  isotropic 

35  36 

materials  have  been  given  oy  Rollins  3  and  Rollins  et  al.  Dunham  and 
37 

Huntington  used  this  technique  for  fused  silica  and  sinqle  crystals 
of  NaCl .  A  treatment  of  the  nonlinear  interaction  of  ultrasonic  waves 
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in  the  framework  of  quantum  theory  has  been  aiven  by  Bajak  and 
Breazeale. ^ 

The  method  of  determining  the  TOE  constants  presented  in  this 
technical  report,  viz.  the  ultrasonic  harmonic  generation  technique, 
involves  a  measurement  of  the  waveform  distortion  of  an  ultrasonic  wave 
as  it  propaaates  throuqh  a  solid.  Due  to  the  anharmonic  behavior  of  the 

solid,  higher  harmonics  are  generated  during  the  passage  of  an  initially 

1  ?.  39 

sinusoidal  ultrasonic  wave  throuch  the  solid.  Breazeale  and  ford 
studied  these  properties,  and  correlated  them  with  the  solution  of  the 

nonlinear  equation  for  a  longitudinal  wave  propaqating  throuah  a  solid. 

40 

Gauster  and  Breazeale  developed  a  capacitive  detector  which  is  capable 
of  absolute  measurement  of  the  amplitudes  of  the  fundamental  and  harmonics 
of  the  waveform  impinging  on  the  end  of  the  sample.  Three  linear  com¬ 
binations  of  the  TOE  constants  of  copper  measured  at  room  temperature 

41 

are  reported  by  them.  The  harmonic  generation  technique  has  the 

advantage  that  it  can  be  applied  to  soft  metals  as  well  as  to  solids 

which  undergo  lattice  slips  and  cracks  on  application  of  pressure. 

Probably  the  greatest  advantage  of  the  harmonic  generation  method 

is  that  it  can  be  used  to  determine  the  TOE  constants  as  a  function  of 

temperature.  The  temperature  dependence  of  the  TOE  constants  is  very 

important  in  the  study  of  anharmonicity  of  solids.  Hydrostatic  and 

uniaxial  pressure  techniques  seem  to  have  practical  problems  at  low 

temperatures.  Harmonic  generation  technique  has  been  used  by  Mackey 
4? 

and  Arnold  to  measure  TOE  constants  of  strontium  titanate  and  by  Meeks 
and  Arnold4^  to  measure  the  temperature  dependence  of  combinations  of 
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TOE  constants  down  to  :  100  K.  Peters  et  al .  developed  a  variable 

qap  capacitive  detector  capable  of  makinq  measurements  at  low  temperatures. 
45 

These  authors  measured  combinations  of  TOE  constants  of  copper  down  to 

liquid  nitrogen  temperatures.  Subsequently  the  apparatus  has  been 
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refined  and  measurements  have  been  made  on  germanium,  fused  silica, 

48 

and  copper  down  to  liquid  helium  temperatures. 

The  work  presented  in  this  report  is  the  result  of  the 
measurement  of  the  ultrasonic  nonlinearity  parameters  and  some  combina¬ 
tions  of  TOE  constants  of  silicon  as  a  function  of  temperature  between 
room  temperature  and  liquid  helium  temperature.  The  apparatus,  experi¬ 
mental  procedure  and  results  are  qiven  in  the  following  chapters.  Our 
measured  nonlinearity  parameters  combined  with  a  well-established  lattice 
dynamical  model  for  diamond-like  solids  has  enabled  us  to  isolate  all 
the  six  independent  TOE  constants  of  silicon  and  germanium  and  plot  them 
as  a  function  of  temperature.  This  is  the  first  report  of  the  temperature 
variation  of  all  the  individual  TOE  constants  of  any  material  as  a  function 
of  temperature.  This  work  is  explained  in  detail  in  the  last  chanter. 

Also  we  have  computed  some  anharmonic  parameters  as  a  function  of 
temperature. 

B.  THEORY 

In  this  section  we  give  the  theory  on  which  the  ultrasonic  wave 
distortion  measurements  are  based.  The  general  equation  of  motion  for 
plane  elastic  waves  propaqatinq  through  a  medium  is  derived,  ft  is  then 
specialized  to  wave  propagation  along  the  symmetry  directions  in  cubic 
crystals.  Pure  mode  longitudinal  wave  propagation  is  possible  along 
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the  [100],  [110]  and  [111]  symmetry  directions  of  a  cubic  crystal. 
Including  the  nonlinear  terms,  one  finds  that  oure  transverse  modes  do 

not  exist.  The  transverse  wave  is  always  accompanied  by  a  lonoitudinal 

50 

wave.  However,  pure  longitudinal  modes  continue  to  exist  for  all 
three  directions  [100],  [110],  and  [111]  even  when  the  nonlinear  terms 
are  included. 

The  wavelength  of  the  sound  waves  in  the  materials  used  in  this 

investigation  are  quite  larqe  compared  to  the  interatomic  spacing. 

Typically  the  wavelength  of  the  ultrasonic  wave  is  about  six  orders  of 

magnitude  greater  than  the  interatomic  spacing.  The  number  of  atoms 

affected  by  one  period  of  vibration  is  of  the  oroer  of  10^.  So  the 

solid  can  be  regarded  as  an  elastic  continuum  and  the  theory  of  finite 
1  ? 

deformations  '  in  the  Lagrangian  formulation  can  readily  be  applied. 

The  eouations  of  motion  governing  finite  deformation  in  cubic  crystal 

8  22 

have  been  given  by  Seeger  and  Ruck,  Bateman  et  al.  and  by  Holt  and 

51 

ford.  The  equations  of  motion  are  derived  here  from  Lagranoe's 
equations,  the  approach  used  by  Holt  and  Ford,  and  then  they  are  aoplied 
pure  mode  longitudinal  wave  propagation  in  cubic  crystals.  From  the 
solution  to  the  resulting  nonlinear  wave  equations  we  derive  expressions 
for  the  ultrasonic  nonlinearity  parameters. 

1 .  The  Equations  of  Motion  for  Plane  Elastic  Waves 
P r onagati ng  through  a  Medium 

In  arriving  at  the  equations  of  motion  of  an  elastic  wave 

propagating  through  a  solid  crystalline  medium,  the  solid  is  considered 

51 

to  be  a  lossless,  homogeneous  and  perfectly  elastic  medium.  The 
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elastic  continuum  is  described  with  the  formalism  given  bv  Murnaahan 
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for  finite  deformations  of  an  elastic  solid.  The  eouations  of  motion 
for  lonq  wavelength  finite  amplitude  sound  waves  in  a  solid  are  derived 
usina  Laqranqe's  eouations  for  continuous  media.  Let  a^ ,  i  =  1,2,3 
be  the  cartesian  coordinates  of  a  point  in  the  unstrained  solid  and  let 
x.j,  i  =  1,2,3  be  the  coordinates  at  time  t  of  the  same  point  in  the 
deformed  solid,  depends  on  the  initial  coordinates  a ■  and  time  t. 

As  has  been  aiven  by  Murnaahan,  let  the  Laqrannian  strain  matrix  n  be 
niven  by 

n  =  \  (J*  J  -  I)  (1.5) 

where  J  is  the  Jacobian  matrix  with  matrix  elements 


JkP  =  3xk/3\  •  (K6) 

J+  is  the  transpose  of  J  and  I  is  the  unit  matrix.  If  p  is  the 
unstrained  mass  density  and  4>(n)  is  the  elastic  potential  energy  per  unit 
unstrained  volume,  then  the  Lagrangian  of  the  system  is 

3 

L  =  i  j  p  X,2  -  4>(n)  .  (1.7) 

c  i=l  1 


Laoranne's  equations  take  the  form 


_d_ 

dt 


d 

dak 


(1.8) 


(1.9) 
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from  Eq.  (1 .5)  wo  have 


’v. 


,  ax.  ,  ax. 

1  ,  _ i  ,  1  ,  _j_ 

2  jk  aa  2  \k  3a, 

^  J 


(1.10) 


■•k  being  the  kronecker  delta.  Substituting  (1.10)  in  (1.9)  we  have  the 
equations  of  motion  given  by 


x . 

i 


3 

=  1 

■  k=l 


_d 

da, 


3 

(  I 

e=l 


Jie 


3(j> 


3n 


kf 


(1.11) 


where 


_2f_  _  _ 3£ 

['r,k?  ’  'n?k 


(1.12) 


This  way  of  derivinq  the  equation  of  motion  given  by  (1.11) 
requires  only  the  definition  of  strain  and  the  assumption  that  the 
elastic  energy  is  a  function  of  strain  alone.  This  eauation  of  motion 
is  derived  for  a  lossless  elastic  continuum  and  is  valid  for  a  real  solid 
under  the  conditions  of  negligible  attenuation  and  dispersion. 

The  equation  of  motion  qiven  by  Eo.  (1.11)  is  to  be  simplified 
to  make  experimental  measurements  and  hence  determine  the  elastic  con¬ 
stants.  Let  us  consider  plane  waves  pronanating  alonq  a  sinqle 
direction  and  orient  the  a-j  axis  along  the  propagation  direction.  Then 
Eq.  (1.11)  simp! if ies  to 


c 


d  y  j  ( n ) 

da,  k£,  ik  3nlk 


(1.13) 


The  elastic  constants  are  defined  as  coefficients  in  an  exnansion  of  the 
elastic  energy  in  powers  of  the  strain  when  the  strain  is  calculated  with 
the  a^  axes  parallel  to  the  symmetry  axes  of  the  crystal.  Thus,  in  order 


w 
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that  On)  in  Eq.  (1.13)  may  contain  the  usual  elastic  constants,  we 
must  determine  the  change  in  $  and  n  under  a  rotation  of  the  ak  axes. 

Let  a  bar  superscript  denote  a  Quantity  calculated  with  the  ak  axes 
parallel  to  the  crystal  symmetry  axes.  Quantities  without  a  bar  super¬ 
script  are  calculated  with  the  a^  axis  parallel  to  the  propagation 
direction. 

If  R  is  the  matrix  which  rotates  the  a^  axis  into  the  a-j  axis 

12 

and  R*  the  transpose  of  R,  Murnaghan  has  shown  that 


n  =  R*  n  R  . 


(1.10 


The  elastic  potential  energy  is  a  scalar  under  rotation,  and  so 


On)  =  <j>(n)  =  OR*  n  R)  , 


(1.15) 


i .e. ,  if 


-,-x  C11  2  .  Clll  3 

On)  ■  y  nn  +  -j-  nn  + 


(1.16) 


then 

On)  =  \  Cn  (.£  Rij*njkRkl  ^  +  F  Clll  (  /k  Rlj*n.]kRkl  ^  +  ■” 

'  (1.17) 

Thus  Eq.  (1.13)  will  contain  the  usual  elastic  constants  if  we  write  it 
in  the  form 


••  _  d  ,  r  ,  3d)(R*  r,  R 

i  da1  k=]  ik  3nlk 


(1.10 


being  the  strain  calculated  taking  the  axis  along  the  nropaaation 


di recti  on. 


If  are  the  components  of  the  particle  displacement, 
Uk=xk-ak,  (1.19) 


the  x^  axes  and  a^  axes  being  parallel.  Substituting  for  <j> ( n )  from 
(1.17)  for  the  elastic  enerq.y  appropriate  to  the  crystal  symmetry  in  (1.18) 
and  performing  the  operations  indicated,  we  obtain 


li,tt  AijUj,aa  +  Bi j*Uj ,aU£, 


aa 


(1.20) 


where 


and 


Uj,a  *  8V3al 
ui,tt  =  ^V8^ 


U.,aa  *  8W 


(1.21) 


In  the  above  equation,  A^-  are  linear  combinations  of  SOE  constants  and 

B...  are  linear  combinations  of  SOE  and  TOE  constants.  Since  we  are 
1  J  K 

interested  only  in  the  lowest  order  of  nonlinearity,  the  terms  involving 

fourth-  and  higher-order  elastic  constants  are  neglected  in  (1.20). 

Let  us  introduce  new  dependent  variables  in  order  to  remove 

the  linear  coupling  of  the  three  compounds  Uk  through  the  terms 

}  A^U.  ,=  •  The  coefficients  A.,  in  (1.20)  form  a  symmetric  matrix  so 
i  j  j  ,  a  a  i  j 

that  an  orthogonal  matrix  S  exists  such  that  S*A  S  is  diagonal.  Hence  if 


we  introduce  the  transformation 


A 


into  Eq.  (1.20),  we  obtain 


p . 

J.tt 


v  -P. 
j  J,aa 


3 

l 


v>.  P„  P 
urn  p  .a  m.aa 
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(1.23) 


where  u •  are  the  eigenvalues  of  the  A..  matrix.  The  v.  form  linear 
J  ij  Jf.m 

combinations  of  SOE  and  TOE  constants.  Equations  (1.23)  can  be 
simplified  further  using  perturbation  theory.  A  straightforward 
perturbation  calculation  shows  that  only  those  terms  in  Eq.  (1.23)  for 
which  j  =  t.  =  m  or  for  which  i  =  m  and  u-  =  u  yield  a  second  harmonic 

J 

which  increases  linearly  with  propagation  distance.  Such  terms  are  said 
to  be  resonant.  Then  the  experimentally  relevant  equations  of  motion 
for  a  plane  wave  take  the  form 


P  •  —  11  .  P  •  + 

j,tt  j  j ,aa 


v...P.  P. 

J  J  J  0,a  j ,aa 


(1.24) 


where  j  =  1,2,3.  A  perturbation  solution  to  (1.24)  will  be  of  the  form 


Bj  s1n<-t  -  kja)  -  (vjjj/8uj>Bj2kj2a 

x  cos(2<i)t  -  2k^a) 


(1.25) 


where  u.  is  the  angular  frequency  and  k.  are  the  wave  numbers  which  are 
related  by 


m 2  =  (vj/o )kj2  .  (1.26) 

a  is  the  propagation  distance  and  B.  is  the  amplitude  of  P.  at  a  =  0. 

J  J 

The  second  harmonic  amplitude  of  Pj  is  given  by 


H2j  =  -<"ijj/8“j>Bj2kjJa  ' 


(1.27) 
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Holt  and  Ford5^  have  tabulated  the  coefficients  and  for  different 

propagation  and  polarization  directions  for  a  number  of  crystals 
belonging  to  cubic  symmetry. 

2 .  Ultrasonic  Nonlinearity  Parameters 

The  ultrasonic  nonlinearity  parameter  is  defined  as  the  negative 
of  the  ratio  of  the  coefficients  of  the  nonlinear  term  to  the  linear 
term  in  the  nonlinear  wave  equation.  As  can  be  seen  from  Eq.  (1.24), 
the  nonlinearity  parameter  can  be  written  as 


=  - ( V  .  .  ./u  ■  ) 
v  JJJ  J 


(1.28) 


The  magnitude  of  the  parameter  8  is  a  measure  of  the  extent  to  which  the 
waveform  becomes  distorted  in  a  specific  propagation  distance,  and  so  is 
a  direct  measure  of  the  nonlinearity  of  the  medium.  It  is  proportional 
to  the  ratio  of  the  second  harmonic  amplitude  to  the  square  of  the 
fundamental  amplitude.  The  parameter  varies  for  different  directions  of 
propagation.  The  expressions  for  these  parameters  for  longitudinal  wave 
propagation  along  symmetry  directions  in  cubic  crystals  result  from 
considering  the  wave  propagation  along  those  directions  as  shown  in  the 
following  sections. 


3 . _ Wave  Propagation  in  Cubic  Crystals 

The  description  of  the  finite  deformation  of  an  elastic  solid 
differs  in  two  major  respects  from  that  of  infinitesimal  deformation. 
Due  to  the  large  deformations  involved,  the  initial  coordinates  of  a 


particle  in  the  undeformed  state  are  not  interchangeable  with  the  final 


coordinates  in  the  deformed  state.  Moreover,  the  expression  for  the 


strain  energy  in  terms  of  the  strains  must  be  changed.  Finite  elastic 
strain  can  be  treated  in  two  different  formulations.  In  the  Lagrangiar. 
formulation,  the  strain  is  described  in  the  undeformed  state  and  the 


initial  coordinates  of  a  material  particle  a^(a,b,c)  are  taken  as 
independent  variables.  In  the  Eulerian  formulation  the  strain  is 
described  in  the  deformed  state  and  the  instantaneous  coordinates  of 
a  material  particle  x^(x,y,z)  are  taken  as  independent  variables.  In 
the  following  we  consider  plane  wave  propagation  along  the  symmetry 
axes  of  a  cubic  crystal  in  the  Lagrangian  formulation. 

a.  Plane  wave  propagation  along  [100]  direction.  Let  us 
consider  first  the  case  of  plane  finite  amplitude  waves  propagating  alona 
the  [100]  direction  of  a  cubic  crystal.  The  displacement  components  in 
this  case  become 


J  = 


a,t)  ,  v  =  v(a,t) 

and 

w  = 

matrix  given  by  (1 

6), 

namely 

1  1 

Tx 

3x 

3x~ 

11  J12  1 3 

3a 

3  b 

3C 

21  J22  J23 

aa 

3b 

3C 

1  i 

3 z 

3Z 

3Z 

31  °32  J33_ 

_3a 

3b 

3C 

— J 

(1.29) 


1  +  u_ 


L w 


ub 

1  +  vt 

Wu 


1  +  w_ 


(1.30) 


where  u 


•  u/'ia,  etc.,  takes  the  form 


52 


J  = 


+  u  0  0 

a 

v  1  0 

a 

w  0  1 

a 
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(1.31) 


x,  y,  z  in  (1.30)  are  the  same  as  Xp  x2,  x3,  respectively,  in  (1.6) 
and  similarly  a,  b,  c  in  (1.30)  are  the  same  as  a-j,  a2,  a3,  respectively, 
in  (1.6).  The  transpose  of  the  Jacobian  is 

J*  = 


1  +  “a  va 


(1.32) 


Substituting  (1.31)  and  (1.32)  into  (1.5)  we  find  that  the  only 
nonvanishing  strain  components  are 

"11  -  “a  +  1  '“a"  *  *a2  *  “a2> 

n12  =  n21  =  \  va  •  ^ ' 

n13  =  n31  =  \  wa 

Substituting  (1.34)  in  the  expression  for  elastic  strain  energy  (1.3) 
and  using  Birch's  values  for  the  TOE  constants,  we  have 


*  =  2  C11 nll 2  +  C44^nl 2  +  n212  +  n31 ^  +  n132^  +  C1 11 nll 3 


+  2  C166^nll ^n122  +  n21  +  n31  +  n132)- 


(1.35) 


differentiating  (1.35)  with  respect  to  the  strains. 
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3<f> 


3n 


11 


Cllnll  +  3Cmn112  +  \  C166(n122  +  n212  +  r,31 2  +  n132] 


3(j> 


3n 


21 


2C44n21  +  C166nlln21 


d<p 


3n 


31 


2C44n31  +  C1 66nl 1 n31 


(1.36) 


Let  us  introduce  the  enqineerinq  stress  tensor  which  is  not  symmetric  as 


T  =  J(^) 

dr) 


(1.37) 


where 


3<t> 

3n 


(1.38) 


Then  the  equations  of  motion  in  the  Laqranqian  coordinate  system  can  be 
written  as 


3Ti  j/9aj  =  Vi  (1,39) 

where  T..  is  the  stress  matrix,  a.  are  the  Laqranqian  coordinates,  t  „ 

1  J  vJ 

is  the  density  in  the  unstrained  state  and  u^  =  xi  -  ai  is  the 
displacement.  The  stress  matrix  can  be  written  in  the  component  form 


as 
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f" 

T1 2 

T1 3 

’ 

"Jll 

J12 

J13 

T21 

T22 

T23 

J21 

J22 

J23 

T31 

T32 

T33 

1 

J31 

J32 

J33 

-J 

J14 

3n 


n 

3<*> 


3ri 


21 


•’■I' 


an 


12 


94?  3  4 


9nj-|  3n 


a* 

3 1) 


3  (J> 

3n00  3n 


13 

34, 


22  "23 


3$ 


32  Sn33 


(1.38) 


from  which  we  can  write  the  stress  components  as 


T„  =  J,,  +  J,, 

11  11  3n^  12 


T12  =  J11  3n^2  +  J12 

Ti 3  =  Jn  +  Ji2 


3$ 


3n 


-  +  J 


3<J» 


21 


a<t> 

3n22 

<}<$> 


13  3n 


+  J 


31 

34, 


13  3n 


3n 


+  J 


32 

dip 


*<t>  +  1 

21  '  "21  3m ,  J22 


T21  =  J 
Tno  =  J 


23 

a<t> 


1 3  3n 


33 


'll 


3o 


+  J 


21 


23  3n 


31 


■7  +  1  9  4  . 

22  "  21  3n,o  J22  ~  J 


34 


'12 


T23  =  J21  3^3  +  J22 

T  =  1  ^  +  1 

'31  J31  •»n11  J32 

T  -  1  _ _ 1  1 

'32  ‘  J31  3 r i -j 2  32 

j  -  1  +  1 

33  J31  Sr,13  J32 


3n22 

** 

:,r|23 

36 


23  3n 


32 


*  ^  r1-  ! 

23  9„33 


+  J 


34 


71 


33  3 


+  J 


'31 

3A 


22 


3ri 


+  J 


33  3r,32  | 

_J*±_  ! 


23 


33  3 n 


33 


(1.39) 


Substituting  from  (1.31)  and  (1.36)  in  Eqs.  (1.39),  we  obtain  the 
following  expressions  for  and  T31 . 


?H 


'll 


Cll[ua  ’1"/^  <#a2  *  "a2)]  4  3c1Hua?  4  l  c166<\>?  4  "a?>' 


f21  C44va  +  Clluava  +  2  C166uava 


r31  =  C44wa  +  Clluawa  +  2  C166uawa 


(1.40) 


where  we  have  neglected  terms  higher  than  quadratic  in  the  displacement 
gradients . 

Since  we  are  considering  the  special  case  of  wave  propagation  in 
the  [100]  or  a  =  a^  direction,  we  require  only  the  three  stress  components 
T^,  T^-j  and  .  In  this  case  the  three  component  eauations  of  motion 
according  to  Eq.  (1.39)  are 


1T11/Ia  =  pou 

1 T  2 1 / 4  a  =  pqV  (1.41) 

^  T  3 1  /  ' 9  =  ‘'  ow 

because  the  stress  components  are  functions  of  -  a  onlv  and  riot  a,, 
or  -  c.  Ui  fferentiatinq  (1.4'])  and  substituting  into  (1.41)  and 
rearranqing  the  terms  we  obtain  the  equations  of  motion  usinq  Birch's 
TOE  constants  as 


0V 

0W 


C 1 1  uaa  (3fll  +  6<lll)ual,aa  +  ((11  +  2  f'16C)(vavaa  +  V'aa 
C44vaa  =  (C11  +  2  C166)(uavaa  +  Vaa} 

C44waa  =  (fll  +  \  C166)(uawaa  +  wauaa) 


(1.4?) 


IJsino  the  relations  between  Birch  and  Bruqner  elastic  coi  .tants  (Iqs.  1.4! 


we  can  write  these  equations  in  terms  of  Bruqger  elastic  constants  as 


P0U  *  C11uaa  (3C11  +  Cm)uauaa  +  ^Cll  +  C16S)(Vaa  +  wawaa5 
PnV  -  C -.v  *=  (C,,  +  C,,fi)(u ,v  +  vu„) 


’0¥  44  aa  '  "166MUavaa  vauaa‘ 

’Ow  '  C44waa  =  ^C11  +  C166)(uawaa  +  wauaa‘ 


?9 


0.43) 


If  we  write 

a  =  C 
P  =  C 


11 


and 


44 


5  -  3C-|  ]  +  C  ni 


y  r  Cin  +  Cn 


'11  1 66 
Eqs .  (1.43)  become 


(1.44) 


V  - 

,uaa 

’  ,uauaa  *  Y<Vaa 

1 

:  > 

O 

llVaa 

=  Y^uavaa  *  Vaa> 

'  ow  " 

uWaa 

=  i(“a"aa  +  Vaa1 

(1.45) 


Since  we  attempt  to  propaaate  a  pure  mode  longitudinal  wave  only. 


we  have 

u  =  u(a,t),  v  =  w  =  0  . 

Then  the  three  equations  of  motion  in  (1.45)  reduce  to 


i  Au  -  *u  =  u  u 
0  33  3  33 


(1.46) 


(1.47) 


Hence  we  see  that  a  pure  mode  lonqitudinal  wave  may  propaaate  in  a 
nonlinear  medium  along  the  [ 1 00 J  direction.  However,  such  a  wave  cannot 
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propagate  without  distortion  and  the  generation  of  hiqher  harmonics. 

Let  us  assume  that  6  <<  a  and  apply  a  perturbation  solution  to 
(1.47)  of  the  form 

u  =  u  "  +  u*  (1.48) 


where  u'  <<  u^  and  use  trial  solutions 
u"  =  A  sin(ka  -  wt)  for  6=0 
u'  =  B  a  sin  2(ka  -  wt)  +  C  a  cos  2(ka  -  wt)  . 
Substituting  (1.48)  in  (1.47),  we  have 


p0u 


mU 


aa 


-  ;ltU 


aa 


(u°  + 


u')a(u°  + 


') 

'aa 


(1.49) 

(1.50) 


(1.51) 


From  the  linear  elastic  equation  of  motion,  we  know  that 

v 0  •  <’-52> 

1  /2 

where  the  phase  velocity  Cq  =  («/oq)  .  Therefore  Eq.  (1.51)  reduces 

to 


»•*’  '  «m’  ■  *  W>  ♦  ua'“aa']  •<’•»> 


Since  u'  u'  we  may  neglect  the  second  and  third  terms  on  the  right 
hand  side  of  (1.53)  in  comparison  with  the  first  term  and  hence  to  a 
first  approximation  the  eauation  of  motion  reduces  to 


OqU  -  ™u 


aa 


Sua  uaa 


(1.5a) 


Substituting  the  trial  solutions  given  in  Eqs.  (1.49)  and  (1.50)  into 
Eq.  (1.54)  we  obtain  the  relation 


?  2  I 

-4pgiii'Ba  sin  2(ka  -  wt)  -  4pQw  Ca  cos  2(ka  -  -nt ) 

-  2»kB  cos  2(ka  -  wt)  -  2ntkB  cos  2(ka  -  wt)  i 

p  1 

+  4>tk£Ba  sin  ?(ka  -  wt )  +  2akC  sin  2(ka  -  wt)  .  (1.55) 

p  ! 

+  2a kC  sin  2(ka  -  wt)  +  4ak  Ca  cos  2(ka  -  wt)  j 

i 

3  , 

=  -5  K-  A2  sin  2(ka  -  wt)  . 

4  I 

Equating  the  coefficients  of  sin  2(ka  -  wt)  and  cos  2(ka  -  wt)  we  find 
that 

-4p  w2Ba  +  4ak2Ba  +  4«kC  -  -(6kV/2)  1 

I  0-56) 

-4pgw2Ca  -  4akR  +  4ak2Ca  =  0  ' 

2 

Recalling  that  a  =  PgCg  from  Eq.  (1.52)  and  that  the  angular  frequency 
w,  the  wave  number  k,  and  the  phase  velocity  Cg  are  related  by  w  =  Cgk, 
Eqs.  (1.56)  reduce  to 

B  =  0,  C  =  -[(kA)27RpgCn2]  .  (1.57) 

Thus  after  one  iteration  we  have  the  approximation  solution 

u(a,t)  =  A  sin(ka  -  wt)  -  [(kA)2<s/8pgCg2]a  cos  2(ka  -  wt)  (1,58) 

which  involves  the  second  harmonic.  If  additional  iterations  are 
performed,  higher  harmonic  terms  will  be  obtained  (but  higher-order 
elastic  constants  would  be  introduced  into  the  expression  for  t(n),  of 
course) . 
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Breazeale  and  Ford  have  applied  the  results  of  finite  amplitude 
wave  analysis,  as  worked  out  for  fluids,  to  pure  mode  lonqitudinal  non¬ 
linear  elastic  waves  in  cubic  sinqle  crystals.  They  write  the  eauation 
of  motion  equivalent  to  Eq.  (1.47)  in  the  form 


1  Outt  =  K2(uaa  +  3uauaa)  +  K3uauaa 


where  for  the  [100]  direction 


(1.59) 


f<2  —  C  -j  -j  ar,d  k3  C 1 1 1 


(1.60) 


using  Brugqer's  notation  for  the  TOE  constants  and  subscript  t  denoting 
time  differentiation.  Comparison  shows  that  Eq,  (1.47)  is  identical 
with  Eq.  (1 .59)  if  we  set 


a  =  K2  and  «  =  3K2  +  K3  .  (1.61) 

2 

If  we  recall  that  oQCg  =  a  from  Eq.  (1.52)  and  substitute  this 
solution  along  with  Eqs.  (1.61)  into  Eq.  (1.58)  we  have 


(3K  +  K  )  ? 

u(a,t)  =  A  sin(ka  -  mt)  -  - ^ -  (kA)  a  cos  2(ka  -  „t)  (1.62) 


which  is  also  the  expression  given  by  Breazeale  and  Ford.  The  nonlinearity 
3  Ko  +  K3 

parameter  e  =  -  - r? - -  can  be  recognized  in  the  second  harmonic 

*2 

amplitude  in  Eq.  (1.62). 

The  equations  of  motion  and  their  solutions  obtained  above  can 


readily  be  specialized  to  wave  propagation  in  isotropic  solids.  Thus 
we  see  that  a  pure  mode  lonqitudinal  nonlinear  elastic  wave  may  nropaqate 
along  the  [100]  axis  or  in  an  isotropic  solid,  but  such  a  wave  generates 
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high'’  harmonics  and  the  waveform  distorts.  A  cure  mode  transverse 
nonlinear  elastic  wave  cannot  propagate  without  the  simultaneous  existence 
of  a  longitudinal  wave. 

b.  Plane  wave  propagation  along  [110]  direction.  Let  us 
consider  the  propagation  of  plane  finite  amplitude  waves  along  the  [110] 
direction,  and  rotate  our  coordinate  axes  from  a  =  [100],  b  =  [010], 
c  =  [001]  to  a'  =  [110],  b'  =  [T10],  c1  =  c  =  [001].  Thus  the  coordinate 
transformation  is  given  by 

a 

-  (R)  b  (1.63)  J 

c 

where  the  transformation  matrix  is 
1//2  1//2  0 

R  =  -1//2  1//2  0  .  (1.64) 

0  0  1 

With  respect  to  the  new  coordinate  axes,  the  displacement  components  are 
u '  =  u ' (a ' ,t)! 

v'  =  v'(a',t)  .  (1.65) 

w'  =  w' (a1 ,t) 

Hence  J'  and  J'*  are  given  by  (1.31)  and  (1.32)  with  the  appropriate 
primes  inserted,  and  the  nonvanishing  Lagranqian  strain  components  in 
tne  primed  coordinate  system  are  of  the  same  form  as  (1.34),  namelv 
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Since  the  expression  for  the  strain  energy  as  given  by  (1.3)  is  written 
in  terms  of  the  unprimed  coordinate  system,  we  must  transform  the  strain 
components  into  the  primed  system  in  order  to  use  equations  (1.66). 
Making  use  of  the  transformation  matrix  (1.64),  we  transform  the  strain 
components  by 


(n)  =  (R*)(n')(R) 


(1.67) 


to  obtain 

nl  1  =  I  ^11 '  '  n12 '  ‘  n21  ‘ ) 

nl 2  =  1  ^nll '  +  nl 2 '  "  n21 '  * 

nl 3  =  1/v^  n7 3 ' 

n21  =  1  ^nll '  '  n12 '  +  n21 ' ^ 

r'22  =  ^  n 1 2 '  *  ^21  ^  (1 .68) 

n23  =  1/v^  nl 3 ' 

n31  =  n31  ’ 

n32  =  n31  ' 

n33  =  0  • 


Substituting  (1.68)  into  (1.3)  we  obtain 
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^  =  4  Cn^n11  '  +  ^ nl 2 *  +  r|21  '  ^  ^ 

+  l  C12^nn  '2  “  ^ nl 2 '  +  n21 '  ^ 

+  1  C44[rin'2  +  (n12'  ■  n2l')2  +  2(n13l2  +  n31 ' 2) ] 

+  l  Cm^nu'2  +  3nl  1 '  ^n12 '  +  n2l'^ 

+  i  C112^nH  '3  ■  nl  T  '  ^n12  ‘  +  n21  '  ^ 

,1_p  ii  i?  ,  i2\ 

+  4  Ci 44  nH  (nl3  n31  ' 

+  \  C1 66^nl  1  ' 3  +  nll '  ^n12 '  ‘  n21  '  )2  +  nll '  (n13‘2  +  n31  '2^ 

+  1  C456n1 1 '  nl 3 ' n31 '  •  (1 -69^ 


Differentiating  (1.69)  with  respect  to  strains,  and  using  the  relations 
nij  =  njf  and  n^-'  =  n  ^  1  we  get 


r)({) 


3nn' 


2  ^Cll  +  C12  +  2CW  nll '  +  4  Clll(nll'2  + 


4n^  ,2) 


+  \  ci 1 2^3n1 1 ’ 2  "  4nl 2 ' 2^  +  2  C144n13 
+  4  C166^3nll  '2  +  2nl 3 '  ^  +  \  C456nl 3 


,2 

.2 


9n 


21 


(C1 1  "  C12^n12 '  +  3C1 1 1 nll ' nl 2 '  "  C1 1 2nl 1 ' nl 2 ' 


—  Op  I  i  1  p  I  |  |  1  n 

3^  ^L44n13  2  L144nl  1  n13  2  L166n11  f,13 


+  lc 


?  u456nn  nl  3  • 


(1.70) 


Substituting  from  Eqs.  (1.31)  with  primes  inserted,  and  from  Egs.  (1.70) 
and  (1.66)  into  Egs.  (1.39)  with  primes  inserted,  we  obtain  the 
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following  expressions  for  T i  ,  T^-j '  and  ' : 


V  =1(C11  +  C12  +  2C44)<uA*  +lu;'2  +  lva'2  +  lw;'2) 

+  I  Cni(ua’2  +  va‘2^  +  4  C112^3ua'  "  va '  ^ 

+  l  C144wa'2  +  I  C1 66^3ua '  +  \  wa'  *  +  8  C456wa' 

T21 '  =  1  ^Cll  +  C12  +  2C44^ua,va’  +  1  ^Cll  "  C12^va' 


+  2  (3Cni  -  C112)u^,w^t 


r3l'  "2  (C11  +  C1 2  +  2C44)ua ,wa '  +  C44wa' 
+  l  (ci44  +  C1 66  +  C456^ua'wa' 


(1.71) 


where  we  have  neglected  terms  greater  than  ouadratic  in  the  displacement 
gradients.  Again  we  use  the  reduced  forms  of  the  component  equations  of 
motion  given  by  Eqs.  (1.41)  with  primes  inserted  since  we  are  only 
concerned  with  waves  propagating  along  a'  =  LUO].  After 
differentiation  of  Fqs.  (1.71),  substitution  into  (1.41)  primed,  and 
rearrangement,  we  obtain  the  equations  of  motion 


P0U  -  aUaa  =  6UaUaa  +  YVaVag  +  Y*WaWaa  | 

V  -  xVaa  =  Y(uavaa  +  Vaa5 
°0^  ‘  x*waa  =  Y*(uawaa  +  wauaa} 


(1.72) 


where  we  have  dropped  the  primes  with  the  understanding  that  in 

Eqs.  (1.72)  u  and  a  are  along  [ 1 1 0 J ,  v  and  b  are  along  [110]  and  w  and  c 

are  along  [001].  Again  we  have  collected  the  linear  terms  on  the  left 


of  the  equality  and  the  nonlinear  on  the  right.  Moreover,  in 
notation, 

n  =  2  ^11  +  ^12  +  2^44^ 

4  =  |  C(c]i  +  ci2  +  2C44^  +  'Clll  +  C1 1 2  +  C166^ 

Y  =  \  ^Cll  +  C12  +  2C44^  +  ^3C1 1 1  "  C112^ 

Y*  =  4  C2 (C-j  1  +  C12  +  2C44)  +  (C144  +  C166  +  C456)J 

x  =  \  (Cll  *  C12^ 

x*  -  c44 

or  in  Brugger's  notation 

4=|  (C17  +  C12  +  2C44)  +  4-  (C171  +  3C-J12  +  12cU6^ 

Y  =  1  ^Cll  +  C12  +  2C44^  +  1  ^Cm  ‘  C11 2 ^ 

y*  =  j  (Cn  +  C12  +  2C44  +  C]44  +  C]66  +  2C456)  . 

Since  we  attempt  to  propagate  a  pure  mode  longitudinal 
Eqs.  (1.72)  reduce  to  the  form  of  (i.47)  qiven  by 

P0U  -  auaa  =  dUgUaa  . 

If  we  rearrange  Eg.  (1.75)  and  use  Eqs.  (1.61)  we  can  rewrite 
the  form  used  by  Breazeale  and  Ford,  namely 

V  ’  *  3Vaa>  *  K3Vaa 
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Birch 1  s 


(1.73) 

(1.74) 

wave, 

(1.75) 
(1.75)  in 

(1.76) 


where  for  the  [110]  direction 
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K2  "  01  2  ^11  +  ^12  +  ^44) 

(1.77) 

«3  =  6  -  3a  =  ^  (Cn]  +  3C]]2  +  1 2C? 66 ) 

using  Brugger's  notation  for  TOE  constants.  Thus,  we  see  that  in  this 
case  also  a  pure  mode  longitudinal  nonlinear  wave  may  propagate  with  the 
generation  of  a  second  harmonic  aiven  by  (1.62),  with  K?  and  substi¬ 
tuted  from  (1.77).  In  addition,  the  waveform  will  be  distorted.  As 
in  the  linear  case,  two  uniquely  polarized  pure  mode  transverse  waves 
can  propagate  but  in  the  nonlinear  case  each  of  these  pure  transverse 
waves  will  be  accompanied  by  a  different  longitudinal  wave. 

c.  Plane  wave  propagation  along  [111]  direction.  Let  us 
consider  the  propagation  of  plane  finite  amplitude  waves  along  the  [111] 
direction,  and  rotate  the  coordinate  axes  from  a  =  [100],  b  =  L010], 
c  =  [001]  to  a'  =  L 1 11 ] ,  b'  =  [112],  c1  =  L110].  Thus,  the  coordinate 
transformation  is  given  by 


a' 

r.i 

b' 

=  (R) 

b 

c* 

c 

) 

where  the  transformation  matrix  is 


/ 2 

Jl 

/2 

R  =  — 
/6 

1 

1 

-2 

-/3 

/3 

0 

With  respect  to  the  new  coordinate  axes,  the  displacement  components 


are 


u 1  =  u ' (a ' ,t) 
v'  =  v' (a1 ,t) 
w '  =  w' (a ' ,t) 
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(1.80) 


Hence  the  nonvanishing  Lagrangian  strain  components  in  the  primed 
coordinate  system  are  the  same  as  those  in  Eqs.  (1.66),  namely 


"11 
nl  2 ' 
nl  3 ' 


"  V  +  2  '-a 


y21 

r31 


l  K.2  +  v:,2  +  w:,2)  ] 

i 


2  V 


1  2  wa 1 


(1.81) 


Since  the  expression  for  the  strain  energy  as  given  by  (1.3)  is  written 
in  terms  of  the  unprimed  coordinate  system,  we  must  transform  the  strain 
components  into  primed  system  in  order  to  use  Eqs.  (1.81).  Making  use 
of  the  transformation  matrix,  Eq.  (1.79),  and  transforming  according  to 
(1.67)  we  obtain 


"11 

1 

6 

(/2 

"21 '  " 

/6 

"13'  + 

1/2  "12* 

/6 

"13'  + 

2"n 

) 

"12 

1 

"  6 

(/2 

"21  ‘  " 

/6 

"31  '  + 

/2  n-j  2  '  + 

/6 

"13 '  + 

2,1  n 

) 

"13 

1 

6 

(/2 

"21  ’ 

-  /6  n31  ' 

2/2  n12' 

+  2"n') 

"21 

1 

6 

(/2 

"21 '  + 

/6 

"31  '  + 

1/2  "12'  ' 

/6 

"13'  + 

2"n 

) 

"22 

1 

'  6 

(/2 

"21 '  + 

/6 

"31  '  + 

/2  n-|  2  1  + 

/6 

"l  3 '  + 

2"l  1 

) 

"23 

1 

6 

(/2 

"21  ’  + 

/6 

"31  '  " 

2/2  n12‘ 

2n 

n') 

"31 

1 

~  6 

(-2/2  n21’ 

*  *  ■'U 

-  S  "13' 

*  2",r> 

1 

r  r  _2_  1 3  ,  1/  ,2  1 2 

C 4 5 6 ^ 2 7  r’ll  9  "11  ”n31  "n 


13  '31  '13 


1  i/  i2  >2  ■  > \ 

9  "11  '*n21  '  "12  "21  "12  ' 


,  {2  ,  o  ,3  ,  ,3 

54  ~2n21  ’  2"l2 


3 " 2 ^  ‘ » i -j  p  '  ^  +  3n i  p  '  -j  '  2  ] 


■^2  ,  i/  ,  2  «  i2  «  i  iNj.  1 1  1 2 

18  "21  ' "l 3  "  2n31  2n31  "13  '  "12  ("31 


2r,l  3  ’  +  2n31  ’  nl  3  ’ 

(1.83) 


Differentiating  Eq.  (1.83)  with  respect  to  strains  and  using  the  symmetry 
relations  n  - i  n  ^  and  n . . '  =  n^-',  we  have 

3^'  *1  (Cll  +  2C12  +  4C4^"ll  '  +  3  Clll  ("11  ?  +  4"21  +  4n31  '2) 

.  1  r  |2  ,  1  r  /  ,2  ,2  ,2> 

3  C112"ll  3  C1 44 ( " 1 1  '  "21  "  "31  > 

3  C1 66  "l  1  "21  "31  ^ 


+  9  ^123("ll  ”  2"?i  "  2"^i  )  +  q 


31  '  9  456v  1 1 


sr 


3  (Cll  "  C1 2  +  2C44)n2l'  +  3  C11 1 (2nl 1 '"21 ' 


r~  1 2  .  py  1 2  ■,  2/2  r  ,  ,2  , 2 , 

/2  r,21  /2  n31  )  +  3  ci  1 2'n21  "  "31  ' 


+  J  C144 (~2nl  1  '"21  '  '  1/2  "21 '  +  1/2  "31 


+  6  ^1 66 (^"l  1  ' n21  '  +  ^  '  _  ^  "ti  '  ) 


+  9  ci23(‘"n  '"21  ^  "21  ’  +  ^  "31 '  ^ 


+  18  C456('2"ll '"21  +  1/2  n21  '  '  ^  n31  '  * 


4 


>  n 


31 


2  (r  _  r  +  r  ) ,  1  +  4  r  t .  > .  '  + 

3  ^  11  1 12  L44;  31  3  L11 1 '  11  31  ‘  21 


4/2 

3 

1  , 


C1 1 2  n  2 1 

1  f 

"31 

144^”"l  1 

1  1 

"31 

+  "2  '21 

"31 '  ^  +  3  1 

166(?'  11 ''  31  '  ‘ 

,  .  '  M 

21  31  ' 

123('"ll 

'  -  • 

'31 

+  2/2  n21 

'"3i ' ' 

456("nll 

1  1 

n31 

+  .2  n21  ' ' 

31  '  -  2/2 

1  ?r  3i ■ 

n  .84) 

Substitutinq  from  Eqs.  (1.31)  with  primes  inserted  and  from  Eqs.  ( 1 . 84 )  and 
from  Eqs.  (1.81)  into  Eos.  (1.39)  with  primes  inserted,  we  obtain 

T1 1 '  =  1  (Cll  +  2C12  +  4C44)(ua'  +  \  ua,?  +  \  va'?  +  2  wa,?) 

1  "  (u' . 2  +  v',2  +  w',2)  +  -  '  2 


+  3  Clllvua 


1 


3  C1 1 2ua  ’ 


,  •  p  /  i  2  1  ,2  1  i  2  \  2 

3  C144(V  "  4  va'  '  4  wa'  *  * 


.  1  r  I  .  2  1.2  1  ,  2.  1  _  /0  ,  2  1  , 

9  C123(V  2  V  ‘  2  V  )  9  C456(?V  '  4  V 


2  r  /  ,  ?  x  1  ,21  ,2 

oElcc(u,  +  .  v  ,  +  .  w  , 

3  166  a  4  a  4  a 

1 

1  a 


o  1  ? 

c  l  |  4  \ 

W  ,  ) 


T21 '  =  3  ^Cll  "  C1 2  +  C44^va '  +  3  (C11  +  2C12  +  aC44 )ua  * 

.  2  r  ,  ,  ,  /2  ,  2  /?  ,  2n 

3  Clll^ua,v  a'  '  4  va '  4  wa '  ^ 


a  1  va ' 


x  /2  r  ,  ,  2  ,  2W  1  _  ,  ,  ,  /?  ,  2  /?  ,  ?s 
6  C112(va'  "  wa'  5  6  Cl44('ua,va'  '  4  V  +  4  wa ’ 


p  /oi  i  ,  ^  2  »  2  *  2  ,  2  v 

1  Va  1  +  Z"  V3  1  *  4  W=,  '  ) 


1 

6  "166vtua'V  '  T  'a 


.  l  r  ,  ,  ,  /2  ,  2  x  /2  ,  ; 

9  C123('ua' V  "  2  va '  2  wa ’ 


1  r  I  ,  i  ,  /2  i  2  »  2  ,  2  \ 

18  C456('W  T  Va '  "  T  wa '  ^ 
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T31 '  =  3  (C11  C12  +  f  44)wa'  +  3  (C11  +  2C12  +  4C44^ua'wa' 


+  3  C1U(2VV  +  /?  VV 


3  C112va'wa' 


+  J  C144('ua,wa'  +  2  va,wa,)  +  3  C166(ua'wa'  "  l  va,wa') 


+  9  ri?3('ua  ,wa'  +  >?  va,wa'^  +  18  C456(-Ua'wa'  "  2~va,wa')  0.B5) 


where  we  have  neglected  terms  higher  than  the  Quadratic  in  the  displacement 
gradients. 

We  can  use  the  reduced  forms  of  the  component  equations  given  by  (1.41) 
with  primes  inserted  because  we  are  only  concerned  with  waves  propaoatino 
along  the  a'  -  LIU]  direction.  After  differentiation  of  Eq.  (1.85), 
substitution  into  Eqs.  (1.41)  primed,  and  rearranqement ,  we  obtain  the 
equations  of  motion  as  follows-. 


0U  ' 

lUaa 

Vaa  +  '(Vaa 

+ 

w  w  )  ' 

a  aa '  < 

ov  ' 

vaa 

'(Vaa  +  Vaa  } 

+ 

1 

f  (w  w  -  V  V  )  1 

a  aa  a  aa  | 

] 

(1.86) 

0W  ‘ 

■  w 

aa 

' (uawaa  +  wauaa} 

+ 

•  (v  W  +  W  V  )  j 

a  aa  a  aa'  1 

where  we  have  dropped  the  primes  with  the  understandi no  that  in  Eos.  (1.86)  u 
and  a  are  along  [111],  v  and  b  are  alona  [112]  and  w  and  c  are  along  LUO]. 

As  before,  we  have  collected  the  linear  terms  on  the  left  of  the  equality  and 
the  nonlinear  on  the  rinht.  In  Birch's  notation  of  TOE  constants 


44 


a  3  ^Cn  +  2C12  +  4C44' 


4  q  [9(C-|-j  +  2  C  i  2  +  40^ )  +  6  C  i  -j  -j  +  120-,-,^,  +  6C 


44 ; 


'in 


'112  1 44 


+  +  2C19,  +  4Cy,rC] 


1  r 


166  1 23  ’  ‘Tl'456- 

111  ”  J,'l  44  '  '•'v'166  "  ‘-"123 


Y  18  ^^11  +  2^12  +  ^^44 )  +  1 2  C  i  i  -  3C1/t/1  +  6C-,cc  -  2C 

1 


C456^ 


(C 


C1 2  +  C44 ' 


f  '  36  (12C111  ’  1 2C1 1 2  +  3C 1 44  3C1 66  +  4C1 23  '  C456) 


(1.87) 


or  in  Brugger's  notation 


6  =  (C11  +  2C1 2  +  4(W  +  9  (C111  +  6C112  +  1 2 C 1 44  +  24C1 66 


+  2C123  +  16C456; 


Y  O  (Cll  +  2C,?  +  4Cy..)  +  Q  (C 


3  '‘'ll  T  tul 2  T  ‘“'44 1  "  9  "'ll!  "  3C1 44  +  6C1 66  "  C123 


2C456) 


Jl 


T-  ~  18  (Cm  "  3 C 1 1 2  +  3C1 44  "  3C1 66  +  2C1 23 


2C456> 


(1.88) 


If  we  attempt  to  propagate  a  pure  mode  longitudinal  wave  only, 
Eqs.  (1.86)  reduce  to  the  form  of  Eo.  (1.47)  given  by 


n0u  ‘  "uaa  =  *uauaa 


(1.89) 


If  we  rearrange  Eq.  (1.89)  and  use  Eos.  (1.61),  we  may  rewrite  Eg.  (1.86) 
in  the  form  used  by  Breazeale  and  Ford,  namely. 


■0“  =  Muaa  +  3u»0  +  u=u 


a  aa 


\  u  u 

3  a  aa 


(1.90) 


I 


where  for  the  [111]  direction 
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^2  ~  n  ~  3  ^11  +  ^12  +  ^44^  I 

K3  =  6  - 3nc  =  ^  (C]11  +  6C112  +  12C144  +  24C166  +  2C123  +  1 6C456) | 

(1.91) 


using  Brugger's  notation  for  TOE  constants.  Thus,  we  see  that  in  this  case 
also  a  pure  mode  longitudinal  nonlinear  wave  may  propagate  with  the 
generation  of  a  second  harmonic  given  by  Eo.  (1.62)  with  and  K3  substituted 
from  (1.91).  Due  to  the  presence  of  e  terms,  transverse  wave  propagation  is 
more  complicated  than  the  previous  [100]  and  [110]  cases. 

Summarizing  the  results  obtained  in  the  above  three  sections  we  can  say 
that  in  the  principal  directions  in  cubic  crystals  the  equations  of  motion 
for  plane  finite  amplitude  longitudinal  waves  are  described  by  equations  of 
motion  of  the  form  (Eqs.  1.47,  1.75  and  1.89) 


n0u 


uU,,  =  fU  U  , 
33  3  33 


or  in  the  form  (Eqs.  1.59,  1.76  and  1.90) 


(1.92) 


n0u  =  K2(uaa  +  3uauaa)  +  K3uauaa 


(1  .93) 


given  bv  Breazeale  and  Ford.  In  the  above  two  equations  the  constants  are 
related  by 


k2  =  ■<  1 

I 

and  |  (1.94) 

i 

K3  -  '•  -  3 1  , 


The  values  of  and  «3  along  the  three  pure  mode  directions  are  given  in 
Table  I  (Egs.  1.61,  1.77  and  1.91). 


46 


Table  1.  Values  of  the  quantities  K2  and  K3  for  longitudinal  wave 

propagation  along  the  pure  mode  directions  in  a  cubic  crystal 


Direction  j 
of  Wave  j 

Propagation  j _ 

[ioo]  ic. 


[110]  (C11+C12  +  2C44)|I  (C111+3C112  +  12C166) 

!  j 

[HI]  j 3  (cn  +C12  +  4C44):1  (Cm  +  6C112  +  12C144  +  24C166  +  ?C123  +  16C456] 


The  ultrasonic  nonlinearity  parameter,  as  defined  in  Section  2,  is  the 
negative  of  the  ratio  of  the  nonlinear  term  to  the  linear  term  in  the  wave 
equation  (1.92  or  1.93),  namely. 


h  =  - 


3K2  +  K3 


(1 .951 


For  an  initially  sinusoidal  disturbance  at  a  =  0,  the  solutions  are  of 
the  form  (Eqs .  1.62) 


A,  2k2ap. 

u  =  A  -j  sin(ka  -  ut)  +  - g—  cos  2(ka  -  cut)  (1.96) 

A, 2k2ap 

where  A,  is  the  fundamental  wave  amplitude  and  A9  =  = - is  the 

C  O 

second  harmonic  amplitude.  A2  contains  the  nonlinearity  parameter  f. 
m  terms  of  A-j  and  A2<  B  is  given  by 


f  '  5  <“2>  S  =  '  J  <3  *  4’  - 


(1  .97) 
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Obviously  a  is  the  propagation  distance  and  k  =  ?n/\ 
So  a  measurement  of  A-|  and  A,,  leads  to  p  which  can  be 
which  are  combinations  of  TOE  constants. 


s  the  wave  vector, 
used  to  evaluate 


CHAPTER  II 


APPARATUS  AND  MEASUREMENTS 

The  ultrasonic  nonlinearity  parameters  are  determined  by  measurinq 

the  amplitudes  of  both  the  fundamental  component  and  the  generated  second 

harmonic  component  of  a  longitudinal  wave  propagating  along  the  three 

symmetry  directions  of  the  sample.  In  the  theory  of  finite  amplitude 

distortion  in  solids  the  ratio  ^/A-j  ),  as  given  in  (1,97),  is  independent 

of  A-j  only  in  the  limit  of  infinitesimal  amplitude.  So  the  smallest 

amplitude  consistent  with  useful  signal  to  noise  ratio  are  used.  The 

proportionality  between  the  second  harmonic  amplitude  and  souare  of  the 

3Q 

fundamental  amplitude  holds  strictly  true  '  for  infinite  discontinuity 
distance  which  is  equivalent  to  infinitesimal  amplitude.  So  measurements 
with  the  smallest  fundamental  amplitude  most  exactly  satisfy  the 
assumptions  made  in  the  theory. 

The  frequency  of  the  fundamental  signal  is  chosen  as  30  MHz  as 
a  compromise  between  the  two  following  considerations:  The  generation 
of  the  second  harmonic  amplitude  is  directly  proportional  to  the 
square  of  the  frequency.  So  higher  frequencies  produce  an  improvement 
in  the  signal  to  noise  ratio.  At  the  same  time,  attenuation  and  effects 
of  nonparallel ism  of  the  sample  faces  also  increase  with  frequency. 

30  MHz  frequency  is  found  to  be  most  favorable.  In  the  nonlinearity 
measurements  the  pulse-echo  technique  is  used,  allowing  one  wave  pocket 
to  die  out  before  the  next  one  is  introduced.  The  pulse  repetition  rate 
is  kept  down  to  about  60  Hz  to  minimize  sample  heating.  The  larqest 
ultrasonic  displacement  amplitudes  attained  at  30  M(,z  are  of  the  order 
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of  10  A°.  From  Eq.  (1.96)  it  can  be  estimated  that  the  second  harmonic 
amplitude  will  be  about  1  percent  of  that  of  the  fundamental.  Thus  a 
detector  is  needed  of  sensitivity  sufficient  to  measure  displacement 

-9 

amplitudes  of  the  order  of  10  cm. 

1.  Room  Tempo ra f ■  m-p  Apparatus 

(a)  The  capacitive  Detector.  A  capacitive  detector  has  been 

40  -12 

developed  to  measure  displacement  amplitudes  as  small  as  10  cm 

(;  10~^  A°).  A  simplified  diaqram  of  the  capacitive  detector  assembly 

is  shown  in  Figure  1.  The  assembly  consists  of  a  detecting  electrode 

held  in  place  by  a  fused  silica  optical  flat  so  that  it  is  insulated 

from  the  outer  ground  ring  which  forms  the  ground  of  the  electrical 

system.  The  detector  and  ground  ring  are  made  optically  flat  to  an 

accuracy  of  half  a  wavelength  of  helium  discharge  light.  The  dete  .or 

electrode  is  positioned  at  the  center  of  the  ground  ring.  It  is 

recessed  slightly  with  respect  to  the  surface  of  the  ground  ring.  The 

sample  faces  are  also  made  optically  flat  to  the  same  accuracy.  The 

sample  rests  on  the  outer  ground  rinq.  The  assembly  is  aligned  properly 

so  that  the  sample  face  and  the  electrode  surface  form  a  parallel  plate 

capacitor.  Typical  gap  spacing  between  the  sample  face  and  the  electrode 

is  of  the  order  of  5-10  microns,  resulting  in  a  capacitance  of  100-75  PF 

for  a  detector  diameter  of  0.916  cm.  A  bias  voltage  of  the  order  of  1 00- 

150  volts  is  applied  to  the  electrode,  producing  very  high  electric  field 

(  1000  kV/cm)  in  the  capacitor  gap.  When  a  plane  longitudinal  wave 

impinges  on  the  sample  face,  causing  it  to  vibrate,  the  gap  spacing  is 

changed  correspondingly  and  an  alternating  voltaqe  is  induced  between 

the  electrode  and  ground.  The  exact  alignment  of  the  sample,  electrode 


Figure  1.  Cross  sectional  view  of  the  room  temperature  apparatus 


button  and  qround  ring  is  very  important  because  a  slioht  misalignment 
can  result  in  shorting  between  the  sample  face  and  the  electrode.  For 
perfect  alignment,  the  electrode  is  fixed  at  the  center  of  the  ground 
ring,  and  an  optical  flat  is  used  in  monochromatic  light  to  produce 
fringes  on  the  ring  and  electrode.  Alignment  of  the  frinaes  assures 
parallelism  of  the  ground  ring  and  the  electrode. 

Onto  the  top  of  the  sample  a  piezoelectric  X-cut  quartz  transducer 
of  fundamental  frequency  30  MHz  is  attached  with  stopcock  grease  as  the 
bonding  material.  The  transducers  used  are  of  diameter  1.27  cm.  A 
copper  electrode  separated  (electrically  insulated)  from  the  outer 
portion  of  the  assembly  by  a  teflon  ring  is  pressed  against  the  transducer 
by  a  sprinq.  The  sprinq  provides  electrical  connection  between  the 
electrode  and  the  BNC  connector  through  which  the  electrical  pulse  is 
applied.  The  entire  system  is  aligned  and  held  in  position  by  steel 
pins.  The  outer  parts  of  both  the  detector  and  the  generator  transducer 
assembly  make  good  contact  with  the  sample  surfaces  so  that  electrical 
shielding  is  obtained  which  is  very  essential  for  accurate  measurements. 

A  pnotograph  of  the  assembled  room  temperature  apnaratus  is  shown  in 
Figure  2. 

b. _ Calibration  Procedure.  Calibration  of  the  capacitive  detector 

and  subsequent  absolute  amplitude  measurements  are  done  by  introducing  a 
substitutional  signal  in  place  of  the  acoustical  siqnal.  Provision  is 
made  in  the  experimental  setup  to  connect  a  substitutional  signal 
generator  and  an  RF  voltmeter  to  the  detector  button.  The  assembly  is 
made  in  such  a  way  that  the  capacitive  detector  need  not  he  removed  from 
the  circuit  during  calibration  measurements.  The  equivalent  circuit  for 
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the  detector  is  the  Norton  equivalent.  The  equivalent  circuit  Tor  the 
detector  is  shown  in  Figure  3.  When  the  acoustic  signal  is  turned  off, 
the  imaginary  switch  is  opened,  and  the  switch  Sp  is  closed  by  connect¬ 
ing  the  substitutional  signal  generator  to  the  detector  assembly.  The 
series  combination  of  the  inductance  L,  the  amplifier  input  impedance 
in  parallel  with  the  detector  capacitance  and  stray  capacitance 
form  the  total  impedance  of  the  detector  circuit.  The  substitutional 
signal  current  i  flows  through  the  same  impedance  as  the  acoustical 
signal  current  i  .  The  inductance  L  is  the  inductance  of  the  wire 
leading  from  the  banana  jack  to  the  BNC  connector  (see  Finure  1).  The 
inductance  of  the  connector  between  the  banana  jack  and  the  bottom  of 
the  detector  can  be  neglected  due  to  the  large  diameter  of  this 
connector. 

If  the  end  of  the  sample  vibrates  sinusoidally  at  an  annular 
frequency  <.>  with  amplitude  FA,  where  A  is  the  acoustic  wave  amplitude 
in  the  sample,  then  the  qan  spacinq  chanqes  with  time  accordinn  to  the 
relation 

S  -•  Sn  +  FA  sin  „t  .  (?.  1  ) 

The  factor  ?  enters  in  the  amplitude  because  the  vibration  amplitude 
of  the  sample  is  twice  the  wave  amplitude  inside  the  sample.  This 
results  from  reflection  of  the  wave  at  the  stress-free  surface  of  the 
sample.  Thus,  the  capacitance  of  the  detector,  considered  as  a  parallel 
plate  capacitor,  responds  as 

l/r  -  l/r  (I  +  ^  sin  t)  (?.?) 

u  -,r) 


Circuit  of  Detector 
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where  Cp  is  the  static  capacitance  of  the  receiver  which  is  given  !>v 


where  a  is  the  area  of  cross  section  of  the  electrode  button,  c  is  the 
dielectric  constant  of  the  medium  (=1  in  our  case  because  the  medium 
is  air)  and  Sq  is  the  static  gap  spacing.  If  Vp  is  the  bias  voltaoe 
applied  to  the  detector,  the  equivalent  voltage  V  of  the  detector  can 
easily  be  obtained  as 


where  A  is  the  amplitude  of  the  ultrasonic  wave.  The  derivation  of 
this  expression  with  the  help  of  an  equivalent  circuit  is  given  in 
Appendix  A1 . 

In  the  equivalent  circuit  qiven  in  Figure  3,  (Ip  and  G<-  represent 
the  ultrasonic  and  substitutional  signal  generators.  Since  the  voltage 
V  is  related  to  ip  by  ip  =  VwCp  where  w  is  the  angular  freouency  of  the 
ultrasonic  wave,  we  have 


2AV.  ...Cn 
b  D 


(?- 5) 


The  substitutional  current  generator  G<.  is  adjusted  to  give  the  same 
output  from  the  amplifier  as  with  the  acoustic  signal.  When  this 
connection  is  satisfied 


’D  7S  ' 


(3.6) 


Practically,  voltages  are  much  easier  to  measure  than  currents 


and  so  the  voltage  across  G<-  is  measured  and  the  current  i^  is  calculated 
from  this  voltage  and  a  knowledge  of  the  impedance  through  which  the 
current  flows.  The  resistor  R  does  not  act  as  a  pure  resistance  at  high 
frequencies.  So  the  impedance  of  the  resistor  R  must  be  measured  at  each 
frequency  used  in  the  calibration  measurements.  These  measurements  are 
made  with  a  vector  voltmeter.  The  sample,  detector  assembly  and  bottom 
plate  are  removed  from  the  apparatus  and  a  50i,  precision  terminator  is 
placed  on  the  "signal  out"  BNC  connector  as  shown  in  Figure  4.  Both 
vector  voltmeter  probes  (with  isolator  tips)  are  placed  at  point  1  and 
the  phase  angle  between  the  signals  is  zeroed  and  the  amplitudes  are 
measured.  Probe  A  of  the  voltmeter  is  then  left  at  point  1  while  probe  B 
is  moved  to  point  2.  The  qenerator  is  readjusted  to  obtain  thp  same 
reading  of  the  A  channel  amplitude  as  before,  and  the  amplitude  of  the  B 
channel  and  the  phase  between  the  probes  are  recorded.  This  is  done  at 
each  frequency  used  in  the  calibration  measurements.  The  impedance  7 
of  the  resistor  R  can  be  calculated  from 


Z 


--1  V 


ijC  +  R  +  jV! 


B1 


VB2e 


VB2e 


it 


(2.7) 


where  is  the  amplitude  measured  by  probe  B  at  point  1,  is  thp 
amplitude  measured  by  probe  B  at  point  2,  C  is  the  stray  capacitance 
which  appears  at  point  2  (this  includes  the  orobe  tip  capacitance),  Rj 
is  the  resistance  of  the  precision  terminator  (measured  to  0.01  with  an 
impedance  bridge),  ,  is  the  angular  freouency,  ;  is  the  phase  annle 
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between  and  Vg^  and  j  =  -1.  A  plot  of  the  real  and  imaginary 
values  of  Z  versus  frequency,  around  the  fundamental  and  second 
harmonic  frequencies  are  shown  in  Figure  5. 

The  substitutional  current  i ^  is  then  calculated  by 


V, 


z  +  [jo'(CD+Cs)  +  50+j(i)LJ 


(2.8) 


where  is  the  voltage  across  the  current  generator  G^.  The  resistor  R 
has  a  value  of  approximately  10  Kn;  therefore  Z  is  much  larcier  than  the 
other  impedances  in  the  aoparatus,  so  that  the  other  impedances  do  not 
have  to  be  known  very  accurately  in  order  to  calculate  i<-  accurately. 
Also,  the  large  value  of  Z  gives  a  large  voltage  V<-  which  can  be  easily 
measured. 


c.  The  Experimental  Setup.  The  block  dianram  of  the  experimental 
setup  for  the  room  temperature  nonlinearity  measurements  is  given  in 
Figure  6.  A  stable  RF  oscillator  is  used  to  drive  a  gated  amplifier. 

The  pulses  are  passed  through  a  matching  network  and  a  30  MHz  bandpass 
filter  which  insures  spectral  purity  of  the  ultrasonic  wave  even  if  the 
quartz  transducer  is  driven  off  resonance.  The  transducer  has  been 
well  bonded  to  the  sample  with  nonaq  stopcock  qrease.  The  signal  from 
the  capacitive  detector  is  fed  to  either  the  30  MHz  or  60  MHz  bandpass 
amplifiers  and  the  amplitudes  of  the  signals  are  measured  accurately 
with  the  boxcar  integrator  and  vector  voltmeter.  The  silicon  samples 
not  being  good  conductors,  do  not  sufficiently  shield  the  RT  pulse  at 
the  transducer  from  the  detector  so  that  a  30  MHz  direct  feedthrough 
pulse  is  received  at  the  detector  hv  radiation  before  the  ultrasonic 
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signal  reaches  the  detector.  This  feedthrough  pulse  is  oftentimes  much 
larger  than  the  echo  of  interest  and  although  it  can  be  separated  in  time 
from  the  echo  to  be  measured,  it  can  overload  the  30  MHz  amplifier  and  the 
boxcar  integrator.  The  recovery  time  of  the  30  MHz  amplifier  and  the 
boxcar  integrator  would  be  too  long  to  allow  proper  measurement  of  the 
first  echo  if  they  are  overloaded  by  the  feedthrough  pulse.  Whenever 
the  feedthrough  is  excessive  to  make  measurements  difficult,  an  RF 
gate  is  used  to  detect  the  acoustic  echo  and  a  MOSFFT  gate  is  used  to 
prevent  overloading  of  the  boxcar  integrator.  Descriptions  of  the  RF 
gate  and  MOSFET  gate  along  with  their  circuits  are  given  in  Appendix  A2. 

2.  Room  Temperature  Measurements 

The  30  MHz  fundamental  and  the  60  MHz  second  harmonic  signals 
have  been  measured  by  the  boxcar  integrator  and  then  the  continuous 
wave  30  MHz  substitutional  signal  is  introduced  at  the  detector.  A 
30  MHz  filter  is  used  between  the  generator  buffer  amplifier  and  the 
detector  to  insure  spectral  purity  of  the  signal.  The  signal  amplitude 
of  the  substitutional  signal  is  adjusted  with  attenuators  so  that  the 
readings  on  the  boxcar  integrator  coincides  with  the  readings  correspond¬ 
ing  to  the  acoustical  signal.  The  60  MHz  substitutional  signal  is 
derived  by  doubling  the  30  MHz  signal  with  a  ring  bridge  mixer  and 
filtering  the  output  with  a  60  MHz  bandpass  filter.  The  circuit  and 
description  of  the  freguency  doubler  are  also  given  in  Appendix  A2.  The 
continuous  wave  substitutional  signals  in  both  cases  are  measured 
accurately  with  an  RF  voltmeter.  The  bias  voltacie  applied  to  the 
detector  is  noted  and  the  gap  spacing  of  the  receiver  is  measured  bv 
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measuring  its  capacitance  with  an  impedance  bridge.  From  these  the 
amplitudes  of  the  fundamental  and  second  harmonic  signals  are  determined. 

The  measurements  have  been  repeated  for  different  values  of 
fundamental  signal  amplitudes,  the  measurements  beinq  done  on  [100], 

[110]  and  [111]  silicon  samples. 

3.  Velocity  Measurements 

In  order  to  determine  the  1C,  values  given  in  Table  1  (p.  46)  or 
the  SOE  constants  of  silicon,  measurements  of  velocity  of  lonoitudinal 
waves  along  the  three  symmetry  directions  have  been  made  usinn  the 
pulse  overlap  technique  with  the  capacitive  receiver  described  before 
as  the  detector.  A  typical  experimental  setup  for  pulse  overlap  technique 
is  shown  in  Figure  7.  The  qated  amplifier  produces  a  series  of  pulses 
that  overlap  in  the  sample.  As  the  freouency  of  the  CW  generator  is 
changed  the  echoes  go  through  interference  maxima  and  minima.  The 
number  of  minima  (or  maxima)  are  counted  and  their  frequencies  noted.  A 
typical  interference  pattern  obtained  using  this  technioue  is  given  in 
Figure  8.  Interference  minima  are  preferred  rather  than  maxima  in  order 
to  minimize  envelope  effects  in  shift  in  the  position  of  the  peaks  due 
to  nonconstant  amplitude  response  of  the  electronic  equipment  as  a 
function  of  frequency.  The  velocity  C  is  then  calculated  from 

r  =  Mk  (o  o) 

where  Af  is  the  change  in  frequency  corresponding  to  a  change  in  peak 
number  An  and  l_  is  the  lenqth  of  the  sample.  The  velocities  of  longitu¬ 
dinal  waves  along  the  three  symmetry  directions  are  given  by 


I  i qurp  p,  •  •  •  obtained  with  the 

pulse  overlap  t  f  :  ‘  :  •> 1  rul^e  trains 

which  interfere  '  • 


C[100]  =  Cn  =  K2  [100] 

C(110]  =  \  (Cn  +  C12  +  ?C44)  =  K2  [110] 
crm]  =  l  (Cn  +  2C12  +  4C44)  =  K2  [111] 
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(2.10) 


So  a  measurement  of  velocities  leads  directly  to  the  K 2  values  qiven  in 

Table  1  (p.  46).  The  measured  values  are  compared  with  the  values 
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reported  by  McSkimin  and  are  found  to  agree  very  well  within  experi¬ 
mental  uncertainties. 


4^ _ Samples 

The  samples  used  in  the  measurements  are  sinqle  crystals  of 
silicon  cut  alonq  the  [100],  [110]  and  [111]  symmetry  directions  having 
lengths  2.5171  cm,  2.5222  cm  and  2.5248  cm,  respectively.  The  sample 
ends  are  made  optically  flat  to  less  than  15"  of  arc  by  hand  lapping 
and  polishing.  The  end  faces  are  then  made  electrically  conductive 
by  coating  a  copper  coatinq  of  approximately  1000  A°  thick  onto  them 
by  vacuum  evaporation. 

5_._  _  The  Cryogenic  Apparatus 

A  general  cross  sectional  view  of  the  cryoqenic  system  is  shown 
in  Figure  9.  The  dewar  used  is  a  conventional  glass  tyne  incorporating 
two  vacuum  chambers  with  silver  coated  surfaces.  The  cryoqenic  fluid 
chambers  are  shown  as  V-|  and  V2  in  the  figure.  For  operation  between 
room  temperature  and  77K,  liouid  nitrogen  is  used  as  the  cooling  aqent 
in  both  the  chambers.  Near  room  temperature,  it  is  convenient  to 
operate  with  liquid  in  V2  only  and  have  dry  nitrogen  gas  in  the  usually 
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Figure  9.  Cryostat  used  for  the  measurement  of  the  temperature 
dependence  of  the  nonlinearity  parameters  of  single  crystals. 
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evacuated  space  between  V-j  and  Vp.  To  operate  between  77  and  3K, 
liquid  helium  is  used  as  the  coolinq  aaent,  maintaininq  liquid  helium 
in  and  keepinq  Vp  filled  with  liquid  nitroqen.  It  has  been  observed 
that  by  pumping  on  the  liquid  nitronen  in  the  inner  chamber,  the 
temperature  can  be  brought  down  to  about  65° K  thus  avoidinq  the  use  of 
expensive  liquid  helium  between  77  and  65K. 

The  sample  assembly  is  surrounded  by  two  stainless  steel  cans. 

The  inner  one  surrounds  the  sample  and  control  assemblies,  and  the 
outer  one  which  is  concentric  with  the  inner  one  provides  a  means  of 
control  ling  the  amount  of  heat  transfer  between  the  system  and  the 
coolant.  The  outer  can  is  in  immediate  contact  with  the  coolant.  The 
space  between  the  two  cans  is  evacuated  to  provide  an  insulatino  jacket 
around  the  inner  can.  The  cans  are  well  polished  to  reduce  radiation 
losses.  The  cans  are  supported  by  three  thin-walled  cupro-nickel  tubes, 
each  of  which  serves  multiple  purposes.  Two  of  the  tubes  have  smaller 
thin-walled  cupro-nickel  tubes  inside  them  to  form  a  coaxial  transmission 
line:  one  for  the  siqnal  input  to  the  transducer  and  the  other  for  the 
output  from  the  detector.  The  tube  sizes  are  chosen  to  make  a  50 
transmission  line.  These  two  tubes  are  also  used  as  vacuum  lines.  The 
tube  that  houses  the  transmission  line  connector  to  the  transducer  is 
the  vacuum  line  for  the  space  surrounding  the  sample  (the  inside  space 
of  the  inner  can).  The  tube  that  houses  the  transmission  line  connector 
from  the  detector  is  also  the  vacuum  lininq  to  evacuate  the  snace  surround¬ 
ing  the  detector  button  of  the  capacitive  detector.  Photoaraphs  showing 
the  cryogenic  apparatus  are  given  in  Fiaure  10. 


Fiqure  in.  The  cryoaenic  apparatus,  (a)  The  complete  apparatus 
(b)  apparatus  with  outer  can  removed;  (c)  apparatus  with  both  cans 
removed . 


The  capacitive  detector  used  for  low  temperature  measurements 
is  similar  to  that  used  in  the  room  temperature  apparatus,  with  the 
additional  feature  that  the  gap  spacing  between  the  detector  button 
and  sample  surface  can  be  controlled  pneumatically.  The  cross  sectional 
view  of  the  variable  gap  capacitive  detector  is  shown  in  Fiqure  11. 

The  lapped  ground  ring  against  which  the  sample  rests  is  undercut  to 
make  it  a  flexible  diaphragm  approximately  0.04  cm  thick.  The  detector 
space  is  evacuated  throuqh  the  tube  which  houses  the  transmission  line 
for  the  signal  from  the  detector.  The  gap  spacing  is  controlled 
pneumatically  by  regulating  the  Dressure  in  the  inner  can.  A  cut-away 
view  of  the  sample  assembly  is  shown  in  Fiqure  1?..  The  capacitive 
receiver,  the  inner  can  and  the  outer  can  are  bolted  in  place  with  indium 
0-ring  seals  so  as  to  have  good  vacuum  in  the  detector  space,  inner 
can  and  outer  can.  Dry  nitrogen  or  helium  gas  can  be  admitted  into  the 
space  surrounding  the  sample  which  helps  to  prevent  unwanted  precipita¬ 
tion  in  the  sample  space.  The  gap  spacing  is  adjusted  pneumatically 
by  letting  in  or  pumping  out  nitrogen  or  helium  gas  in  the  sample  space. 
The  variable  gap  capacitive  detector  makes  it  possible  to  have  a  constant 
detector  capacitance  for  the  detector  when  +  e  temperature  is  varied. 

The  effects  due  to  differential  therm-’  ox  ''.  ion  of  the  sample  and  the 
associated  assembly  are  thus  compensated  by  regulating  the  pressure  in 
the  inner  can. 

Once  the  sample  has  been  cooled  to  liquid  nitrooen  or  liquid 


helium  temperature,  the  temperature  can  be  increased  with  an  electrical 
heater.  Two  sets  of  heater  coils  are  wound,  one  above  and  one  below  the 
sample  so  that  uniform  heating  of  the  sample  takes  place.  The  temperature 
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figure  12.  Cut-away  view  of  the  low  temperature  sample  assembly. 


! : 

ot  t  hr  '.(imp  1  o  is  measured  with  a  germanium  or  platinum  sensor.  A 
uermanium  sensor  with  negative  temperature  coefficient  of  resistance 
is  used  below  ?0  K  and  a  platinum  sensor  with  positive  temperature 
•  oef '  i i i ent  of  resistance  is  used  at  higher  temperatures.  The 
e  1  t*i  t  r  i  c:a  1  output  from  the  sensor  controls  the  power  delivered  to  the 
neator,  this  being  accomplished  with  an  Artronix  temperature  controller, 
'hr  temperature  controller  contains  a  bridge  circuit.  One  of  the  arms 
o’  the  hridoe  is  the  temperature  sensor.  If  the  bridqe  is  unbalanced, 
the  controller  provides  more  or  less  power  to  the  heater  which  tends 
to  balance  the  bridge.  The  deviation  of  the  temperature  from  the  set 
point,  and  the  voltage  on  the  heater  are  monitored  on  a  dual  pen  strip 
chart  recorder.  Monitoring  the  heater  voltaqe  allows  one  to  determine 
when  thermal  equilibrium  has  been  established  in  the  sample  chamber. 

The  system  allows  one  to  increase  the  temperature  slowly  and  measure 
it  accurately.  A  block  diagram  of  the  temperature  measurement  and  control 
system  is  given  in  Figure  13. 

To  facilitate  accurate  temperature  measurements,  the  measurements 
are  done  with  different  sensors  than  the  temperature  control  sensors. 

A  constant  current  generator,  shown  in  Tinure  14,  supplies  a  current 
to  the  sensor  (either  platinum  or  germanium),  and  the  voltage -is 
measured  across  the  sensor.  A  four-wire  connection  is  made  to  the 
sensors  so  that  the  voltaqe  measurement  is  made  by  a  path  which  does 
not  carry  current.  This  prevents  lead  resistance  from  affecting  the 
sensor  resistance  measurement.  A  voltmeter  which  has  an  input  impedance 
of  1000  nr.  is  used  for  measui  ement  so  that  neglinible  current  is 
drawn.  The  voltmeter  is  also  used  to  measure  the  current  by  measurine 


Voltmeter 


the  voltaae  across  a  known  resistance  in  series  with  the  current -carrying 
lead.  The  constant  current  generators  is  required  because  the  heater 
current  also  flows  throuqh  the  ground  of  the  apparatus.  When  a  simple 
resistive  network  is  used  to  supply  current,  the  noisy  heater  current  is 
larqe  enough  to  change  the  voltage  with  respect  to  qround  which  chanqes 
the  sensor  current  bv  a  noticeable  amount.  The  resistance  measurement 
can  sull  be  made,  in  principle,  as  the  sensor  voltage  changes  when  the 
sensor  current  is  chanqed,  but  since  the  measurements  cannot  be  made 
simultaneously  the  heater  noise  makes  measurements  difficult.  The 
single  transistor  constant-current  generator  is  very  effective  in 
holding  the  current  constant  regardless  of  voltage  changes  with  respect 
to  qround.  The  voltaqes  measured  are  of  the  order  of  1  mV  so  that 
contact  potentials  and  thermal  potentials  need  be  taken  into  account. 

An  effective  way  to  do  this  is  to  take  two  or  more  resistance  measure¬ 
ments  using  different  currents.  Then  the  sensor  resistance  R  is  found 
by 

v  -  v 

R  (Ml) 

?  “  ‘l 

where  V-j  end  are  two  voltages  and  I-j  and  I ?  are  the  corresponding 
currents . 

6 .  C ry o genic  jM on linear ity  Measurements 

The  block  diagram  of  the  apparatus  used  for  low  temperature 


ultrasonic  nonlinearity  measurements  is  shown  in  figure  15.  The  signal 
path  is  the  same  as  that  for  room  temperature  measurements  with  the 
exceptions  that  (1)  the  substitutional  signal  capability  is  absent 


Block  diaoram  of  the  cryogenic  nonlinearity  measurements. 


// 


because  only  relative  amplitudes  need  be  measured  and  {?)  a  Ml  Mil/ 
broadband  filter  is  introduced  between  the  capacitive  detector  and  the 
60  MHz  amplifier  in  addition  to  themostat  gates.  The  feedthrough 
problem  is  much  worse  in  the  cryogenic  measurements  than  in  the  room 
temperature  measurements.  In  the  room  temperature  apparatus,  the 
sample  could  be  adjusted  manually  until  proper  sample  seating, 
principally  on  the  transducer  end,  is  achieved  resulting  in  the  reduction 
of  feedthrough  signal.  In  the  low  temperature  apparatus,  the  sample  is 
inaccessible  to  such  adjustments  and  so  the  two  gates  are  important  in 
the  cryogenic  measurements  to  insure  that  the  amplifiers  are  not  over¬ 
loaded.  Under  these  severe  conditions  the  6D  MHz  amplifier  can  be 
overloaded  by  the  30  MHz  feedthrough.  So  the  60  MHz  filter  is  used  to 
attenuate  the  feedthrough  while  passing  the  60  MHz  signal.  Since  the 
filter  is  a  completely  passive  network,  harmonics  cannot  be  generated 
by  it.  In  Figure  16  a  photograph  of  the  entire  low  temperature  system  is 
given. 

The  low  temperature  nonlinearity  measurements  are  done  using  a 
comparison  technique  because  the  30  MHz  and  60  MHz  amplifiers  have 
gains  which  can  be  functions  of  the  signal  amplitude.  The  procedure 
is  to  set  the  capacitance  of  the  detector  to  a  chosen  value.  As  the 
temperature  is  varied,  the  variations  in  the  capacitance  (or  spacina) 
is  nullified  very  conveniently  by  adjusting  the  pressure  in  the  inner 
can  with  a  regulator  valve.  The  bias  voltage  on  the  detector  is  set 
to  some  predetermined  value  150  volts).  The  RT  power  applied  to  the 
quartz  transducer  is  adjusted  to  obtain  some  predetermined  reading  of 
the  30  MHz  signal  on  the  boxcar  integrator.  The  detector  bias  voltage 
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is  then  chanqed  until  the  60  MHz  second  harmonic  signal  reaches  the 
predetermined  boxcar  integrator  readinq.  The  pulse  repetition  rate, 
pulse  length  and  qains  of  the  30  MHz  and  60  MHz  amplifiers  should  be 
kept  undisturbed  during  the  entire  experiment.  In  order  to  avoid  the 
trouble  of  considering  the  frequency  denendence  of  the  eouipment, 
the  frequency  of  the  CW  generator  is  also  kept  constant. 

Always  the  bias  voltage  of  the  fundamental  siqnal  remains  the 
same  and  any  deviation  of  the  fundamental  signal  amplitude  from  the 
predetermined  boxcar  integrator  reading  is  compensated  by  adjustinn 
the  input  signal  to  the  transducer.  The  second  harmonic  bias  voltane 
is  then  varied  to  obtain  the  predetermined  second  harmonic  boxcar 
integrator  reading.  So  the  only  variable  in  the  entire  temperature 
cycling  is  the  second  harmonic  bias  voltage.  The  output  of  the  60  MHz 
amplifier  is  proportional  to  the  amplitude  of  the  second  harmonic 
and  the  detector  bias  voltage  V^.  If  this  output  is  held  constant  and 
the  measurements  are  made  at  two  temperatures  T-j  and  T^, 


(?.l?) 


where  the  different  subscripts  indicate  the  temperature  corresponding 
to  the  value.  Therefore 


(3.13) 


The  measurements  are  made  relative  to  the  room  temperature  values.  The 
relative  nonlinearity  parameters  of  silicon  along  the  three  symmetry 


80 

directions  have  been  measured  between  room  temperature  and  3°K  by  this 
technique. 

7.  Transducer  Bonding  for  Low  Temperature 
Nonlinearity  Measurements 

The  nonlinearity  measurements  are  made  with  X-cut  30  MHz  transducers 
of  diameter  1.27  cm  bonded  to  the  sample.  For  room  temperature  measure¬ 
ments  nonaq  stopcock  grease  is  used  as  the  bonding  material.  This  bond¬ 
ing  material  fails  at  low  temperatures.  So  for  low  temperature  measure¬ 
ments  a  cellulose  tape  adhesive  is  used  as  the  bonding  material.  A 
cellulose  tape  manufactured  by  Technical  Trade  Corp.,  Carbondale, 

Illinois,  under  the  designation  7510-551-9818,  LT90C  Type  1,  Class  A  is 
used  for  the  purpose.  The  adhesive  is  applied  as  follows.  A  piece  of 
plastic  film  with  a  hole  slightly  larger  in  diameter  than  the  trans¬ 
ducer  is  centered  over  the  sample  end  and  a  piece  of  the  cellulose  tape 
is  applied  which  sticks  to  the  sample  throuqh  the  hole.  The  tape  is 
smoothed  down  to  remove  all  air  bubbles  and  a  drop  of  water  is  applied 
to  the  back  of  the  taoe.  When  the  tape  turns  cloudy  (3-5  minutes),  the 
water  is  blotted  away  and  the  plastic  film  is  gently  lifted.  The  backing 
of  the  tape  separates  cleanly  from  the  adhesive  everywhere  except  at 
the  edge  of  the  hole.  The  adhesive  stretches  between  the  sample  face 
and  the  plastic  film  at  the  edqe  of  the  hole  and  is  very  careful lv 
severed  with  a  sharp  blade  so  that  the  backing  is  oently  lifted  away 
without  disturbing  the  adhesive  bonded  to  the  sample.  All  traces  of 
moisture  on  the  bond  are  removed  by  keeping  the  sample  in  vacuum  for 
approximately  three  hours.  The  transducer  is  then  carefully  placed 
on  the  adhesive  and  all  air  bubbles  are  removed  by  firm  but,  careful 


rubbing  with  the  eraser  at  the  end  of  a  pencil.  This  bond  is  found  to 
work  very  well  over  the  entire  temperature  range. 


CHAPTER  III 


RESULTS  AND  DISCUSSION 

In  this  chapter  we  report  the  results  of  the  measurements  made 
on  silicon  samples  to  determine  their  ultrasonic  nonlinearity  parameters 
at  room  temperature  as  well  as  the  temperature  dependence  of  the  non¬ 
linearity  parameters  between  room  temperature  and  liquid  helium 
temperature.  Temperature  variation  of  some  more  combinations  of  TOE 
constants  are  also  determined  and  reported.  The  strain  generalized 
Gruneisen  parameters  of  silicon  calculated  alonq  the  symmetry  directions 
from  nonlinearity  parameters  are  also  evaluated  and  presented. 

a . _ Room  T emperature  Nonlinearity  Measurements  of  Silicon 

The  results  of  the  absolute  amplitude  measurements  at  room 
temperature  for  silicon  alonq  the  three  principal  directions  are  qiven 
in  Table  2.  In  order  to  satisfy  the  theoretical  assumption  of  infini¬ 
tesimal  amplitudes,  these  measurements  are  made  at  low  amplitudes  which 

give  good  signal  to  noise  ratio.  Irrespective  of  the  value  of  A-| ,  it 

2 

may  be  noted  from  the  table  that  the  value  of  A^/A^  is  a  constant.  The 
variation  of  the  absolute  amplitude  of  the  second  harmonic  with  the 
square  of  the  fundamental  amplitude  for  all  the  three  orientations  are 

plotted  in  Figure  17.  The  straiqht  lines  show  a  well-defined  linear 

2 

variation  of  A£  with  A^  .  The  lines  do  not  pass  throunh  the  origin 
because  of  the  residual  noise  of  the  experimental  setup.  The  slope  ot 
the  lines,  in  the  least  square  sense,  is  the  best  fit  to  the  value  of 
the  quantity 
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Table  2.  Amplitudes  of  Ultrasonic  Wave  Components  for  Silicon  at  Room 
Temperature 


i  Second  ; 

1 

! 

Sample 

Orientation 

Frequency 

Used 

(MHz) 

Fundamental 

Amplitude 

A]  x  10~10  M 

|  Harmonic 

1  Amplitude 

A0  x  10"13  M 

1  a?/a3 
!  z 

x  107  M'1 I 

■  3  A,2  k 

7.1677 

16.5192 

0.321536 

0.6657 

7.0557 

15.7156 

0,315683  1 

0.6536 

6.7944 

14.7334 

0.319154  1 

0.6608 

6.5704 

!  12.9475 

0.299918 

0.6210 

[100] 

30.41 

6.0105 

12.3224 

0.341094 

1  0.7062 

5.6745 

10.8937 

0.338314 

0.7005 

5.3385 

9.3758 

0.328980 

0.6812 

4.5918 

6.9649 

0.330331 

0.6839 

4.2558 

5.8933 

0.325384 

0.6737 

3.6585 

_ 

4.7325 

0.353577 

0.7321 

i 

3.1279 

6.0989 

0.623370 

1.5303 

3.0614 

5.9663 

0.636579 

1 .5184 

2.9282 

5.3034 

10.618519 

1.5997 

2.7508 

4.9056 

0.648297 

1.5628 

[no] 

30.24 

2.7064 

4.7731 

0.651653  1 

1.5944 

2 . 6842 

4.7731 

0.662477  1 

1.5915 

2.5955 

4.3753 

0.649480  1 

1.6263 

1 

2.3958 

3.5798 

0.623674  1 

1  .5311 

2.1296 

2.8373 

0.625618  ! 

1.5358 

2.1075 

2.7843 

0.626875  ! 

j 

1.5389 

3.4860 

5.9676 

0.4910717' 

1  .2555 

3.4442 

5.8433 

0.492585  ! 

1.2594 

1 

3.3190 

5.3460 

0.485305  1 

1.2408 

3.2563 

5.1719 

0.4877545' 

1.2470 

[111] 

30.30 

3.2150 

5.0973 

0.4931491 1 

1  .2608 

3.0267 

4.4757 

0.488565  ! 

1.2491 

2.9015 

4.1027 

0.48733161 

1.2460 

1 

2.7345 

3.6054 

0.482167  ' 

1.2328 

i 

2.6301 

3.3568 

0.4852672! 

1.2407 

i 

2.3590 

2.7352 

'0.491520  1 

1.2567 

! 

! 

2.1918 

2.3622 

:  0 .491  7165! 

1.2572 

30  40  50  60 

A,(io"20  m2 ) 


Figure  17.  Variation  of  with  A]  at  room  temperature  for 


s  i  l  icon . 


8!) 


3  2 

i  1J'k  a 


9  3K9  +  K, 


(3.1  ) 


which  involves  the  nonlinearity  parameter  p.  The  role  played  by  the 

power  lost  from  the  fundamental  to  the  second  harmonic  is  examined  in 

3  2 

Appendix  3.  If  this  effect  is  taken  into  account,  the  quantity  „  pk  a 
is  given  by 


3 

8 


,  2 

k  ap 


(3.2) 


2 

rather  than  A^/A-j  .  For  the  entire  range  of  amplitudes  qiven  in  Table  2, 

2 

the  difference  between  A^/A^  and  the  quantity  calculated  from  Eq.  (3.2) 
is  negligible. 

In  Figure  18  we  have  plotted  the  dimensionless  parameter  . 
calculated  for  various  values  of  the  fundamental  amplitude.  The  scatter¬ 
ing  of  points  in  Figure  18  is  well  within  experimental  error  in  view  of 

the  fact  that  the  measurements  depend  on  absolute  measurements  of  displaco- 

.3 

ment  amplitudes  of  the  order  of  10  A.  In  Table  3  we  tabulate  the  values 
of  ?,  obtained  by  a  least  squares  fit  to  the  data  from  all  the  plotted 
points  for  all  the  three  principal  directions.  The  standard  deviation 
from  the  least  square  fit  calculated  by 


n^/ 

1  i 


-  1 


(3.3) 


is  included  in  the  :  values. 


87 


Table  3.  Values  of  the  Nonlinearity  Parameter  i3  (Least  Square  fit) 


Sample 

Orientation 

Sample 

Length 

(cm) 

[100] 

2.5171 

7.0336 

[no] 

2.5222 

4.6887 

[HI] 

2.5248 

3.7500 
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The  parameters  Kg,  which  are  combinations  of  TOE  constants,  alona 
the  three  directions  are  evaluated  from  the  B  values  by  Eq.  (1.95) 

Kg  =  (3  +  S)K2  (3.4) 

and  are  tabulated  in  Table  4.  The  Kg  values  used  in  the  calculation 

are  evaluated  from  the  SOE  constants  of  silicon  reported  by  McSkimin  an.H 
54 

Andreatch.  We  have  measured  the  velocities  and  verified  that  our 

values  of  SOE  constants  agree  very  well  with  those  of  McSkimin  and 
54 

Andreatch.  The  Kg  values  used  in  the  evaluation  of  Kg  are  also  given 
in  Table  4.  The  room  temperature  TOE  constants  of  silicon  have  been 
measured  earlier  by  hydrostatic  and  uniaxial  pressure  technique  by 

pc 

McSkimin  et  al.  The  Kg  values  calculated  from  their  TOE  constants 

are  also  tabulated  in  Table  4  for  comparison.  It  may  be  noted  that  the 

Kg  values  obtained  from  these  two  entirely  different  techniques  agree 

55 

very  well  within  experimental  inaccuracies.  Keating  has  theoretically 
evaluated  all  the  six  TOE  constants  of  both  silicon  and  germanium  in 
terms  of  three  anharmonic  first  and  second  neighbor  force  constants  and 
two  harmonic  force  constants  on  the  basis  of  central  potential.  We  have 
calculated  the  Kg  values  from  these  theoretical  TOE  constants  and  these 
values  are  also  presented  in  Table  4.  The  agreement  between  experimental 
and  theoretical  Kg  values  shows  that  Keating's  model  is  a  very  good 
approximation  for  silicon  and  other  diamond-like  solids. 

b.  Strain  Generalized  Gruneisen  Parameters  of  Silicon 

In  this  section  we  derive  the  relationship  between  the  strain 
generalized  Gruneisen  parameters  and  ultrasonic  nonlinearity  parameter 


Table  4.  The  K2  and  K3  Values  of  Silicon  Along  the  Principal  Directions 
at  Room  Temperature 


Sampl e 
Orientation 

K9  x  1011 
(dynes/cm  ) 

|  K3  x  1012  dynes/cm^ 

Present 

Experiment 

Experimental 
Values  of 
McSkimin  et  al . 

Theoretical 
Values  of 
Keatinqb 

[100] 

16.5779 

-  8.34  ±0.11 

-  8.25 

-  8.21 

[no] 

19.4470 

-14.95  ±0.22 

-14.75 

-15.68 

[111] 

20.4031 

-13.77  ±0.12 

-13.31 

-13.58 

Reference  26. 


^Reference  55. 
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for  solids.  Since  the  medium  is  not  stressed  in  the  determination  or 
nonlinearity  parameter  by  harmonic  qeneration  technique,  the  qeneralized 
Grtineisen  tensor  calculated  from  the  nonlinearity  parameters  is  purely 
i sen tropic. 

Let  us  rewrite  Eq.  (1.24)  in  the  form 


a2p  u. 

—J-  =  -1  (1  + 
at2  p0 


vjjj  \  8  PJ 
pj  8al  aa^ 


(3.5) 


When  the  ultrasonic  nonlinearity  parameter  R..  defined  by  (1.28)  is  zero, 
Eq.  (3.5)  reduces  to  a  linear  differential  eouation  describinn 
infinitesimal  amolitude  wave  propagation  of  velocity 

cj =  sV/2  •  (3'6) 

For  finite  amplitude  waves  r.  is  nonzero  and  it  is  convenient  to  define 

23 

the  natural  wave  velocity  by 


W(j,N)  =  o  +  ^  ^-)31/2  *  C-0  -  Pi  •  (3.7) 

Pq  Pj  J  J  ‘'a-J 

So  Rj  is  a  direct  measure  of  the  nonlinearity  of  the  solid  with  respect 
to  elastic  wave  propagation. 

The  isentropic  generalized  Grilneisen  parameter  YpS(j,N)  describing 
the  frequencies  u>(j,fi)  are  defined  by 
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s 

y 

rs 


(j.N)  = 


T3TNT 


a«o(J,w) 


3n 


s,n=0 


rs 


(3.8) 


where  "n  refers  to  the  original  coordinate  system  and  s  is  the  entropy. 
The  strain  derivatives  in  the  original  coordinate  system  may  be 
expressed  in  terms  of  the  derivatives  with  respect  to  the  displacement 
gradients  apj/aa-j  in  the  transformed  coordinates  by  using  chain-rule 
differentiation.  Hence,  we  may  write 


( — — ) 

;s,n=0 

3n 


8npq  3(auk/3al) 


rs 


3n 


3n 


rs 


pq 


(3.9) 


Using  Eqs.  (1.14)  and  (1.22)  and  the  relation 


3(3u,/3a,)  , 

(  3^pq  }s,n=0  =  2  (5kp6lo  +  5kq5lp^ 


(3.10) 


where  6..  are  Kronecker  deltas,  we  can  express  Eq.  (3.9)  in  the  form 

*  J 


^s,n=0  2  ^kj^kr^ls  +  RlrRks^a(ap./aa1 ) ^s,3p ./aa^O  ' 

rs  J  J 


(3.11) 


Substituting  (3.11)  into  (3.8)  we  obtain 


Y^J.N)  - 


2  6k^RkrRls  +  RlrRks^w(j^)  x  3 ( apj/aa1 ) ^s.apj/aa^O' 

(3.12) 


In  the  Debye  model,  the  lattice  vibrational  frequencies  w(j,l5)  are  related 
for  any  state  of  strain  to  the  natural  velocity  by^6 


«  W(j ,N)  . 

Hence,  using  Eqs.  (3.7)  and  (3.13)  we  may  write 
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(3.13) 


/  I  MM)  \  =  t  1  /  3W(j.N)  s  , 

'wCMT  afapj/aa^^s.apj/aa^O  W(j,N)  va(3pyaa1  rs.apyaa^O  =  -  ^  '  • 

J 

(3.14) 

Substituting  (3.14)  into  (3.12)  we  obtain  the  general  relationship 
between  the  (y,s)  component  of  the  isentropic  strain  generalized  Griineisen 
tensor  and  the  solid  nonlinearity  parameter  B  •  to  be 

J 

Vs»-N>  MjVVis  +  RisRks>  ■  <3-,5> 


From  this  equation  we  find  that  the  (1,1)  components  of  the  isentropic 
strain  generalized  Griineisen  tensors  for  the  case  of  longitudinally 
polarized  ( j=l )  waves  propagating  in  the  symmetry  directions57  are 


[100])  =  ^  B-j  for  the  [100]  direction 

-I  ( 1 »  [110])  =  j  for  the  [110]  direction,  and 

y^i  (1 ,  [111])  =  ^  &i  for  the  [11]  direction. 


(3.16) 


In  Table  5  we  list  room  temperature  values  of  the  (1,1)  components 

of  the  experimentally  measured  isentropic  strain  generalized  Griineisen 

tensors  for  the  pure  mode  directions  of  silicon  and  germanium  crystals. 

The  corresponding  mixed  state  strain  generalized  Griineisen  parameters 

58 

parameters  calculated  from  the  theory  of  Brugger  using  experimental 


values  of  elastic  constants  are  also  given  in  the  table.  For  silicon 
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the  TOE  constant  data  of  McSkimin  et  al.  and  our  measured  nonlinearity 

parameters  are  used.  For  germanium  the  TOE  constant  data  of  Bateman 

2.2.  46 
et  al.  and  nonlinearity  parameters  data  of  Bains  et  al.  are  used. 

It  must  be  pointed  out  that  exact  agreement  between  the 

i^entropic  strain  generalized  Gruneisen  parameters  y^s  and  the  Brugger- 

B 

Gruneisen  parameters  y“$isnot  to  be  expected  because  the  former  are  of 
a  pure  thermodynamic  state  whereas  the  latter  are  of  a  mixed  thermo¬ 
dynamic  state.  The  difference  between  the  two  may  be  of  the  same  order 
as  the  difference  between  adiabatic  and  isothermal  elastic  constants. 

An  examination  of  Table  5  reveals  that  for  the  [100]  direction  for  both 

s  B 

silicon  and  germanium  the  agreement  between  y^.s  and  y°s  is  good.  They 
agree  reasonably  well  for  the  [110]  direction  but  do  not  agree  for  the 
[111]  direction  for  either  silicon  or  germanium.  Further  work  in  this 
direction  is  underway  to  trace  the  reason  for  the  disagreement  in  the 
results  along  the  [111]  direction. 

c.  Temperature  Dependence  of  the  Nonlinearity  of  Silicon 

The  measured  values  of  |b|  down  to  approximately  3°K  are  given  in 

Table  6  along  with  the  ratio  of  J 6 ]  to  the  room  temperature  |b|.  The 

room  temperature  value  is  the  most  probable  average  value  obtained  from 

Table  3.  These  data  for  all  the  three  directions  are  listed.  The 

values  calculated  from  e  are  also  tabulated  in  Table  6  between  300  and 

3°K.  In  Figure  19  we  have  plotted  the  data  obtained  for  the  parameter 

for  the  three  orientations  as  a  function  of  temperature.  The  S0E 

constants  of  silicon  as  a  function  of  temperature  have  been  plotted  by 
53 

McSkimin  from  room  temperature  to  77°K.  We  measured  the  velocities 


Temperature 

T°K 

k2 

xlO12  dyn/cm2 

Bias  Voltage 
of  the  Second 
Harmonic 

Vb 

i  JLi 

Ip,  1 

I  B  | 

"  K3 

xlO12  dyn/cm2 

[100]  direction 

3  10 

1.6755 

150.01 

1.1086 

0.7509 

8.8005 

4  02 

1.6755 

150.05 

1.1087 

0.7508 

8.8004 

6  00 

1.6755 

150.15 

1.1079 

0.7503 

8.7994 

10  57 

1.6755 

150.85 

1.1028 

0.7468 

8.7803 

16  14 

1.6755 

152.51 

1.0908 

0.7387 

8.7396 

20  11 

1.6755 

153.00 

1.0873 

0.7363 

0.7275 

24  93 

1.6755 

153.00 

1.0873 

0.7363 

8.7275 

35  08 

1.6753 

153.55 

1.0833 

0.7336 

8.7129 

47  82 

1.6753 

155.12 

1.0725 

0.7263 

8.6/62 

59  25 

1.6752 

158.80 

1.0477 

0.7095 

8.5913 

66  33 

1.6750 

160.88 

1.0343 

0.7004 

8. 5445 

77  60 

1.6747 

163.21 

1 .0195 

0.6904 

8.4927 

86  84 

1.6743 

163.45 

1.0177 

0.6892 

8.4847 

101  57 

1 . 6740 

163.58 

1.0163 

0.6882 

8.4781 

119  92 

1.6732 

163.82 

1.0140 

0.6867 

8.4666 

131 .51 

1  6723 

164.08 

1.0118 

0.6852 

8.4545 

147  89 

1.6713 

164.05 

1.0113 

0.6849 

8.4479 

167  99 

1.6700 

164.22 

1 . 0093 

0.6835 

8.4343 

181  11 

1.6689 

164.25 

1 . 0083 

0.6828 

8.4253 

197  93 

1.6674 

164.35 

1.0067 

0.6817 

8.4122 

216.08 

1 .6654 

164.62 

1.0051 

0.6807 

8.3971 

231.91 

1 .6637 

164.41 

1 . 0042 

0.6800 

8.3850 

250.29 

1.6617 

164.42 

1.0031 

0.6793 

8.3715 

266.30 

1 .6601 

164.41 

1.0022 

0.6787 

8.3604 

272.48 

1.6594 

164.45 

1.0017 

0.6784 

8. 3554 

Room  temp. 

1.6566 

164.48 

1 . 0000 

0.6772 

8.3353 

[110]  direction 


2.98 

1.9641 

4.02 

1.9641 

4.90 

1.9641 

8.41 

1.9641 

10.57 

1.9641 

14.00 

1.9641 

16.14 

1.9639 

24.93 

1.9639 

35.08 

1.9637 

47.82 

1.9635 

56.82 

1.9634 

64.58 

1.9633 

77.79 

1.9629 

150.05 

1.2742 

150.08 

1.2740 

152.55 

1.2534 

155.65 

1.2284 

156.48 

1.2219 

157.85 

1.2113 

158.38 

1.2072 

161.22 

1.1858 

165.55 

1.1547 

170.08 

1.1241 

174.66 

1.0946 

177.12 

1 .0796 

178.65 

1.0704 

1.9922 

17.6286 

1.9913 

17.6256 

1.9591 

17.3915 

1.9200 

17.2055 

1 . 9098 

17.1454 

1.8933 

17.0482 

1.8869 

17.0087 

1 .8534 

16.8114 

1.8048 

16.5234 

1.7570 

16.2401 

1.7109 

15.9686 

1.6874 

15.8285 

1 .6730 

15.7405 

I 
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Table  6  (continued) 


Bias  Voltage 


of  the  Second 

_  y 

Temperature 

K2 

12  2 

Harmonic 

1  J_, 

'  ft  ' 

1 6 1 

K3 

12  : 

T°K 

xlO  u  dyn/crri 

vb 

R 

xlO  dyn/cm 

[110]  direction 

(continued) 

84.02 

1.9627 

179.92 

1.0627 

1 .6610 

15.6682 

94.75 

1.9623 

180.89 

1.0564 

1.6512 

15.6073 

107.28 

1.9618 

181.55 

1.0520 

1.6443 

15.5628 

118.77 

1.9614 

182.46 

1.0464 

1.6355 

15.5078 

129.77 

1 . 9606 

182.78 

1.0441 

1.6319 

15.4803 

136.76 

1 . 9602 

183.21 

1.0414 

1.6277 

15.4525 

147.89 

1.9594 

183.56 

1.0389 

1.6238 

15.4232 

160.87 

1.9585 

184.52 

1.0329 

1 .6144 

15.3609 

173.35 

1.9575 

185.21 

1.0284 

1.6074 

15.3120 

184.10 

1.9563 

185.98 

1.0235 

1.5997 

15.2574 

194.91 

1.9552 

186.38 

1.0206 

1.5952 

15.2224 

206.38 

1.9542 

186.78 

1.0179 

1.5910 

15.1900 

223.37 

1.9522 

187.24 

1.0145 

1.5857 

15.1434 

241.70 

1.9501 

187.98 

1 . 0094 

1.5777 

15.0803 

260.13 

1 . 9480 

188.66 

1.0074 

1.5746 

15.0460 

272.48 

1 . 9464 

188.86 

1.0031 

1.5678 

14.9939 

Room  temp. 

1 . 9433 

189.18 

1.0000 

1.5630 

14.9420 

[111]  direction 

3.97 

2.0601 

150.12 

1.0575 

1.3219 

14.3500 

4.90 

2.0601 

150.08 

1 .0577 

1.3221 

14.3513 

7.03 

2.0601 

150.03 

1.0581 

1.3226 

14.3544 

10.08 

2.0601 

149.89 

1.059' 

1.3239 

14.3624 

17.00 

2.0600 

150.26 

1.0565 

1.3206 

14.3413 

24.50 

2.0600 

150.58 

1.0543 

1.3179 

14.3246 

30.70 

2.0599 

150.81 

1.0526 

1.3158 

14.3109 

36.99 

2.0599 

150.98 

1.0513 

1.3141 

14.3004 

43.55 

2.0598 

151.77 

1.0459 

1.3074 

14.2583 

50.47 

2.0596 

152.55 

1.0406 

1 . 3008 

14.2162 

59.25 

2.0595 

152.76 

1.0394 

1.2993 

14.2062 

66.33 

2.0593 

152.81 

1.0390 

1 .2988 

14.2018 

70.15 

2.0592 

153.81 

1.0323 

1.2904 

14.1492 

77.79 

2.0592 

155.00 

1.0245 

1.2806 

14.0886 

86.84 

2.0587 

155.62 

1 .0199 

1.2749 

14.0500 

92.49 

2.0585 

155.82 

1.0185 

1 .2731 

14.0375 

103.85 

2.0581 

156.02 

1.0166 

1.2708 

14.0206 

114.16 

2.0576 

156.44 

1.0136 

1.2670 

13.9937 

129.77 

2.0566 

156.68 

1.0113 

1.2641 

13.9690 

139.68 

2.0561 

156.48 

1.0121 

1.2651 

13.9718 

147.89 

2.0553 

157.08 

1 . 0080 

1.2600 

13.9349 

157.32 

2.0546 

156.83 

1.0093 

1.2616 

13.9400 

167.99 

2.0537 

156.45 

1.0111 

1.2639 

13.9481 

174.54 

2.0532 

156.45 

1 .0109 

1.2636 

13.9429 
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Table  6  (continued) 


Temperature 

T°K 

K2 

12  2 
xlO  dyn/cni 

Bias  Voltage 
of  the  Second 
Harmonic 

Vb 

1  JLi 

br 

M 

'  K3 

12  2 

xlO 'c  dyn/crri 

[HI] 

direction  (continued) 

181.71 

2.0525 

156.45 

1 .0105 

1.2631 

13.9350 

197.93 

2.0508 

156.48 

1.0092 

1.2615 

13.9137 

212.44 

2.0492 

156.48 

1 . 0086 

1.2608 

13.8985 

227.03 

2.0476 

156.75 

1.0061 

1.2576 

13.8680 

244.15 

2.0457 

156.66 

1.0057 

1.2571 

13.8520 

260.13 

2.0436 

156.89 

1 .0034 

1.2543 

13.8207 

272.48 

2.0421 

156.82 

1.0032 

1.2540 

13.8087 

Room  temp. 

2.0388 

157.11 

1.0000 

1.2500 

13.7619 
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and  hence  the  values  at  different  temperatures  belo,.  77  K  and  found 
that  these  values  agree  with  the  curves  of  McSkimin  extrapolated  to  0nK. 
For  calculations  uf  as  a  function  of  temperature  the  values  have 
been  evaluated  as  a  function  of  temperature.  They  are  also  tabulated  in 
Table  6. 

From  Figure  19  it  may  be  noted  that  the  values  undergo 
noticeable  temperature  variation  below  77  K  especially  [110].  Between 
room  temperature  and  3°K,  the  values  for  the  [100],  [110],  and  [111] 
directions  vary  by  5.3?,  15.2?,  and  4.1?4,  respectively.  The  data  for 
the  [100]  and  [ill]  directions  are  quite  smooth  and  are  nearly  parallel 
to  the  temperature  axis.  The  variation  is  maximum  for  the  [110] 
direction  between  77  and  3  K,  the  variation  being  as  high  as  11?.  Since 
B  is  a  direct  measure  of  the  anharmonicity  (nonlinearity),  a  large  change 
in  e  (and  therefore  K^)  at  low  temperatures  is  not  surprising  in  view  of 
the  thermal  expansion  data  which  shows  that  there  is  a  large  change  in 
the  anharmonicity  of  silicon  and  other  diamond-like  solids  below  liquid 
nitrogen  temperature  (thermal  expansion  becomes  negative). 

The  TOE  constant  C-j^  appears  in  the  values  along  all  the 
directions.  It  can  be  subtracted  out  from  [110]  and  [111].  In 
this  way  we  can  get  two  more  combinations  of  TOE  constants,  namely 
(C^2  +  66 ^  anc^  ( C 1 23  +  6^144  +  SC^g)  and  their  temperature  variation 

can  be  studied.  The  temperature  variation  of  these  two  combinations 
along  with  that  of  C-j-ji  are  plotted  in  Figure  20.  As  ran  be  seen  from 
Figure  20,  the  variation  of  these  combinations  as  a  function  of 
temperature  is  not  very  great  down  to  77  K.  But  between  77  and  3  K, 
there  is  a  considerable  variation  in  some  of  them.  In  particular,  the 
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Figure  20.  Variation  of  some  combinations  of  TOE  constants  of 
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combination  ( C i ^3  +  ^^144  +  ®^456 ^  crosses  zer0  at  about  8°K  and  becomes 
positive  below  that  temperature.  Since  these  are  the  TOE  constants  which 
would  play  a  significant  role  in  the  propagation  of  transverse  modes,  the 
lack  of  agreement  between  the  theoretically  calculated  value  of  thermal 
expansion  from  elastic  data  and  measured  thermal  expansion  values  miaht 
be  attributed  to  the  theoretical  assumption  that  the  TOE  constants  are 
independent  of  temperature. 


d.  Temperature  Dependence  of  the  Nonlinearity  of  Germanium 

In  order  to  supplement  the  necessary  data  for  the  work  to  be 

described  in  the  next  chapter,  we  reproduce  here  the  data  on  the 

temperature  variation  of  the  «3  values  of  germanium.  In  Table  7  we 

present  the  measured  values  of  the  nonlinearity  parameter  fl  and  calculated 

values  of  of  germanium  between  300  and  3°K.  These  data  are  adopted 

59 

from  the  Ph.D.  dissertation  of  J.  A.  Bains.  In  Figure  21  we  have 
plotted  the  temperature  variation  of  the  parameters  K3  of  germanium 
along  the  pure  mode  directions  between  room  temperature  and  3°K. 
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Table  7.  Measured  Values  of  of  Germanium  as  a  Function  of 
Temperature  (Ref.  59; 


Sample 

Orientation 

Temperature 

°K 

T2 

K,  x  10 

*  2 

dyn/cm 

2.94 

-  8.47 

9.01 

-  8.47 

19.30 

-  8.47 

27.81 

-  8.47 

43.28 

-  8.46 

56.72 

-  8.46 

68.66 

-  8.46 

77.40 

-  8.22 

100.51 

-  8.29 

[100] 

124.35 

-  8.34 

148.48 

-  8.38 

160.53 

-  8.36 

184.52 

-  8.28 

196.55 

-  8.30 

220.51 

-  8.30 

245.16 

-  8.17 

270.27 

-  7.88 

Room  temp. 

-  7.78 

3.00 

-17.21 

9.01 

-16.32 

19.30 

-15.89 

27.81 

-15.54 

43.28 

-14.91 

56.72 

-14.35 

68.66 

-14.57 

77.68 

-14.62 

[no] 

100.51 

-14.54 

124.34 

-14.45 

148.43 

-14.46 

160.50 

-14.27 

184.57 

-14.24 

196.52 

-14.13 

220.51 

-14.15 

245.11 

-14.20 

270.28 

-14.15 

Room  temp. 

-14.10 

3.01 

-13.28 

9.11 

-13.61 

19.38 

-13.84 

27.78 

-13.72 

43.22 

-13.59 

56.78 

-13.56 

68.69 

-13.28 

103 


Table  7  (continued) 


Sample 

Orientation 


Temperature 

°K 


K,  x  in 

j  2 

dyn/cm 


12 


77.48 

-13.52 

100.53 

-13.31 

124.33 

-13.71 

148.46 

-13.49 

160.53 

-13.40 

184.54 

-13.53 

196.58 

-13.44 

220.51 

-13.47 

245.15 

-13.51 

270.28 

-13.73 

Room  temp. 

-13.40 

[HI] 

(continued) 


dvnes/cm  ' 
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Figure  21.  Temperature  variation  of  the  parameters  of 
germanium. 


CHAPTFR  IV 


r 


TEMPERATURE  VARIATION  OF  THE  TOE  CONSTANTS 
OF  SILICON  AND  GERMANIUM 

The  SOE  constants  of  a  number  of  solids  have  been  measured  as  a 

function  of  temperature  by  a  number  of  authors.  The  first  measurement 

53 

of  this  kind  has  been  done  by  McSkimin  on  silicon,  germanium  and 
fused  silica  from  room  temperature  to  liquid  nitrogen  temperature.  Later 
the  measurement  has  been  extended  to  liquid  helium  temperature  by  Fine^ 
on  oermanium.  Subsequently  measurements  have  been  reported  on  a  number 
of  solids.  The  results  have  been  reviewed  by  Hearmon61  and  are  not 
reproduced  here. 

The  measurement  of  the  TOE  constants  as  a  function  of  temperature 

is  much  more  difficult  than  in  the  case  of  SOE  constants.  If  one  uses 

the  McSkimin  technique,  chanaes  in  sound  velocity  need  be  measured 

with  the  application  of  pressure  and  as  a  function  cf  temperature.  Due 

to  experimental  difficulties  not  much  work  has  been  done  in  this 

29 

direction.  Salama  and  Alers  measured  the  TOE  constants  of  cooper 

CO 

at  295,  77  and  4.2°K.  Sarma  and  Reddy  measured  the  TOF  constants  of 
nickel  at  298  and  80°K.  In  both  cases  the  technique  aave  results  in 
which  the  inaccuracies  of  the  measurement  entered  at  every  temperature 
at  which  measurements  were  made.  The  effect  of  temperature  on  TOF 
constants  are  found  to  be  comparatively  small.  The  normal  tendency  is 
for  the  stiffness  to  increase,  i.e.,  to  become  less  negative  as  the 
temperature  rises.  Some  authors  have  reported  the  results  of  calculations 
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based  on  theoretical  considerations.  These  are  also  reviewed  by 

Hearmon.^  As  already  has  been  pointed  out  in  the  introduction,  the 

ultrasonic  harmonic  generation  measurements  can  readily  be  undertaken 

at  low  temperatures  but  until  now  this  technique  allows  one  to  measure 

only  three  combinations  of  TOE  constants  as  a  function  of  temperature 

and  not  individual  TOE  constants. 

In  this  chapter  we  have  combined  the  results  of  our  measurements 

of  the  nonlinearity  parameters  of  silicon  and  germanium  with  an 

established  lattice  dynamical  model  to  arrive  at  expressions  for 

individual  TOE  constants  in  terms  of  nonlinearity  Darameters  and  have 

studied  their  variation  with  temperature.  The  model  used  is  the  central 
63  55 

potential  model  ’  introduced  by  Keating  for  diamond-like  solids.  For 
the  first  time  in  literature  we  have  studied  the  temperature  dependence 


of  the  individual  TOE  constants  of  a  solid  between  300  and  3°K.  The 
highlights  of  Keating's  theory  are  given  in  the  next  section  and  in  the 
following  section  the  Keating  force  constants  are  expressed  in  terms  of 
the  K^  parameters  along  principal  directions  and  expressions  for  TOE 
constants  are  arrived  at.  Section  c  is  devoted  to  results.  Graphs 
showinq  the  temperature  variation  of  the  TOE  constants  of  silicon  and 
germanium  are  plotted. 

As  is  well  known,  the  diamond-like  solids  silicon  and  nermanium 


exhibit  an  anomalous  thermal  expansion  behavior  at  low 
64,65,66 


temperatures . 


Thermal  expansion  of  silicon  is  negative  between 


120  and  20nK  and  is  again  positive  below  that  temperature.  Thermal 

expansion  of  germanium  is  neqative  between  40  and  16°K  and  is  positive 

68 

below  that  temperature.  The  works  done  have  been  reviewed  by 
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Carr  et  al.  The  Gruneisen  parameter  follows  a  similar  behavior  as 

thermal  expansion.^  Theoretical  calculations  based  on  the  quasiharmonic 
71  72 

approximation  ’  using  room  temperature  values  of  TOE  constants  do  not 
predict  this  behavior. We  have  evaluated  the  Gruneisen  parameter 
at  low  temperatures  using  our  temperature  dependent  TOE  constants  for 
both  silicon  and  germanium.  Results  are  compared  with  previous  publi¬ 
cations.  Temperature  variation  of  some  more  anharmonic  pronerties  like 
pressure  derivatives  of  the  SOE  constants,  the  Anderson-Grunei sen 
parameter, etc.  are  also  studied  and  the  results  are  reported.  The  final 
section  constitutes  discussion  and  some  concluding  remarks. 

a .  Keating's  Model  for  the  Lattice  Dynamics  of  Diamond-Like  Sol  ids 
63  55 

The  Keating  ’  model  for  the  diamond-like  solids  is  basically 

1 3 

equivalent  to  the  Eorn-Huang  ‘  approach  of  imposing  the  invariance 
requirements  on  the  elastic  strain  enerqy  of  the  crystal.  Keatinq's 
theory  provides  additional  insight  into  the  interatomic  forces  in  this 
group  of  solids  and  has  more  operational  advantages  over  the  Born-Huano 
method.  The  method  demonstrates  that  all  purely  first-neighbor  inter¬ 
actions  are  central  onlv.  Keatinq's  approach  confirms  that  in  the 
harmonic  approximation  there  are  no  noncentral  purelv  first-neiahbor 
interactions  present  in  any  nonmetallic  crystal.  In  the  harmonic 
approximation  treatment  to  derive  expressions  for  SOE  constants, 

Keating's  model  involves  only  two  types  of  interactions,  a  nearest- 
neighbor  central  term  and  a  noncentral  second  neighbor  term. 

The  basic  unit  cell  of  the  diamond  structure  is  a  rhombohedron 
with  two  atoms,  atoms  1  and  0  in  Figure  22,  on  its  major  axis  which  is 


and  filled  circles  represen 
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directed  along  the  [111]  direction.  The  three  neighboring  unit  cells  of 

interest  contain  atoms  2  and  5,  3  and  6  and  4  and  7,  respectively.  The 

expression  for  the  microscopic  strain  energy  in  terms  of  the  central 

first  neighbor  constant  «  and  the  noncentral  second  neiohbor  constant  p 
63 

has  been  written  and  after  imposing  crystal  symmetry  reouirements , 
the  expression  in  the  harmonic  approximation  becomes 
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(4.1) 


where  a  is  the  lattice  parameter  and  e..'s  are  the  strain  parameters 

*  J 

with  e^  =  exx  +  e  +  ezZ.  u',  v',  w*  are  the  atomic  displacement  compo¬ 
nents.  By  imposing  the  condition 


au*  3v'  aw'  "  u 


(4.2) 


n  n 


one  obtains 


u* 

=  -a 

V. 

evz 

v’ 

=  -a 

F, 

ezx 

and  w1 

=  -a 

C 

exy 

where 

f.  =  (u-f)/(a+p)  .  ^ ‘ ^ 

After  substitution  for  u\  etc.,  and  comparing  the  resulting  expression 

with  the  well-known  macroscopic  expression  for  strain  energy  for  cubic 
75 

crystals  ,  namely 


u  =  \  cn(exx2  +  eyy2  +  O  +  C12(eyyezz  +  ezzexx  +  exxeyyJ 


+  2  C44^yz  +  ezx  +  ex.y  ^  ’ 


(4.5) 


we  obtain 


Cn  =  (■«  +  3P)/4a 

C-j 2  =  -  (;)/4a  (4‘6) 

^•44  =  (,t P- / a  ( ■  x  +  •>)  ! 

These  form  the  expressions  for  the  SOE  constants  in  terms  of  the  harmonic 
force  constants  a  and  P.  It  can  easily  be  seen  that  the  parameter  i  is 
given  by 


EC.  9 

n  -  P  _  12 

<t  +  p  TCi  1  +  C-|  2 ) 


(4.7) 


A  calculation  of  the  TOE  constants  of  silicon,  oermanium  and 
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other  diamond-like  solids  has  been  presented  by  Keatino  which  is  an 
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extension  of  the  above  method  of  setting  up  the  expression  for  the 
strain  energy  including  the  anharmonic  terms.  Expressions  for  the  six 
TOE  constants  are  derived  in  terms  of  three  anharmonic  first  and  second 
neighbor  force  constants  and  the  two  previously  determined  harmonic 
force  constants.  As  has  been  done  before,  expressions  for  the  strain 
energy  are  set  up  in  the  microscopic  and  macroscopic  reoimes.  Includinq 
the  anharmonic  terms,  the  expression  for  macroscopic  strain  energy  can 
be  written  as 


6  Clll^el  +  e2  +  e3  ^  +  2  C112^el  (e2  +  e3 ^  +  e2  ^e3  +  el  ^ 

+  e32(el  +  e2)]  +  C1 23  ele2e3  +  1  C144(ele4?  +  ^  +  e3e62) 
+  1  C166[el2(e52  +  e62)  +  e2(e62  +  e42)  +  e3<eA2  +  e^)] 

+  C456  e4e5e6  +  \  C1 1 3  +  e2  +  e3  ^  +  2  C12^el  ^e2  +  e3^ 


+  e22(e3  +  e] )  +  e-^e-j  +  e2)l 


(4.8) 


where  e,  =  exx,  e2  =  eyy,  e3  =  ezz,  e4  -  exy,  e5  -  exz,  and  e&  -  eyz> 

Considering  the  microscopic  enerqy  density,  it  should  be  pointed 
out  that  the  most  important  forces  between  atoms  in  diamond-like  solids 
are  apparently  short-range  forces  mainly  due  to  shell-shell  and  shell- 
core  interactions.  Cochran^  has  obtained  a  good  fit  between  theoretical 
and  experimental  phonon  dispersion  and  elasticity  data  with  the  shell 
model  and  only  short  range  forces.  The  shell  model  reduces  to  a  riqid- 
ion  formulation  if  only  the  elasticity  is  to  be  considered  and  excellent 
agreement  has  been  obtained  between  theory  and  experiment  usino  only  a 


nearest  neighbor  interaction  and  a  noncentral  second  neinhbor  interaction . 
Since  the  shell-shell  and  shell-core  effects  give  only  small  contributions 
for  third  and  more  distant  neighbors,  we  include  interactions  out  to 
second  neighbors  only. 

Including  the  anharmonic  terms,  the  expression  for  the  microscopic 

55 

strain  energy  can  be  written  as 

=  [(y1  -  +  3f')/3  +  (a  +  3?)/8aQ](ei^  +  +  e^) 

+  2 ( y '  +  36'  -  fr')e-|e2e3  +  [y '  -  V  +  r'/3  +  (a  -  p)/8ap][e-|2(e2  +  e3) 

+  e22(e3  +  e-j)  +  e3?(e1  +  e2)]  +  2y'(l  -  03e4e5e6 

+  [y'(1  -  r)2  +  f'(  1  +  k)2  +  1  e'(l  +  f.)(3 f  -  1) 

+  -2(i  -  ?)/8aQ][e1e42  +  e2eg2  +  e3eg2) 

+  [Y-(l  -  ,)2  .  o  ’  ( 1  +  r)2  +  1  r  *  ( 1  +  0(3  -  O 

+  ;2(<  -  0/83^^  (e62  +  eg2)  +  e2(e2  +  e4?)  +  e3(e42  +  eg2)]  . 

(4.9) 

•'  and  '  are  the  anharmonic  force  constants.  The  details  of  the 
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derivation  arp  not  given  here  as  it  is  alreadv  available  in  literature.' 
rompar<ng  the  expressions  (4.8)  and  (4.9),  the  expressions  for  the  TOF 
constants  of  diamond-like  solids  are  obtained  as 
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2 

where  y  =  2y  ‘ ,  6  =  24 1  and  t  =  ^  r  ’  are  the  redefined  anharmonic  force 
constants. 

Keating's  model  has  been  used  to  evaluate  the  TOE  constants  of 
silicon  and  germanium  and  it  is  found  that  the  results  are  in  very  oood 
aqreement  with  the  room  temperature  experimental  values.  This  agreement 
is  quite  remarkable  owing  to  the  fact  that  the  model  involves  only  three 
adjustable  parameters  and  it  suogests  that  the  anharmonic  force  constant 
model  used  is  a  fairly  realistic  one.  The  introduction  of  additional 
interactions  into  the  strain  enemy  will,  of  course,  allow  a  much  better 
fit  but  such  a  fit  then  miaht  not  be  sianificant  physically.  The  fact 
that  Keating's  model  involves  only  three  anharmonic  parameters  has 
enabled  us  to  isolate  all  the  TOE  constants  of  silicon  and  qermanium 
from  our  measured  ultrasonic  nonlinearity  parameters  and  study  their 
variation  as  a  function  of  temperature.  The  Keatinq  model  has  been 
applied  by  Nandanpawar  and  Rajagopalan^  to  predict  the  temperature 
dependence  of  C^gg  and  30^  +  4C^gg  between  73  and  ?93°K  of  germanium. 
Their  results  agree  reasonably  well  with  available  experimental  values. 
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b.  TOE  Constants  in  Terms  of  the  Kq  and  Parameters 

The  expressions  for  the  and  parameters  which  are 
combinations  of  SOE  and  TOE  constants  given  in  Table  1  are  reproduced 
here  for  the  ease  of  further  derivations.  These  parameters  are  obtained 
directly  from  harmonic  generation  experiments. 

K^[100]  =  Cn  ] 

.  I 

K2[110]  =  2(C7 1  +  C12  +  2C44)  |  (4.11) 

K2[lll]  =  \(CU  +  2C12  +  4C44)  j 

K3[100]  ~  C-j  •]  •) 

K3(n°]  =  1^111  +  3C1 1 2  +  1 2C1 66^  1 

l 

K3[l 1 1 ]  =  1(C111  +  6C1 1 2  +  1 2C1 44  +  24C166  +  2C1 23  +  16f456^  1 

(4.12) 

Substituting  for  the  SOE  and  TOE  constants  from  Eqs.  (4.6)  and 
(4.10)  in  (4.11)  and  (4.12),  respectively,  we  can  express  and 
along  the  principal  directions  in  terms  of  the  harmonic  force  constants 
a  and  8  and  the  anharmonic  force  constants  y,  4  and  c. 

K2[100]  =  (a  +  3p)/4a 

K2[110]  =  I[(«  +  3p)/4a  +  (a  -  p)/4a  +  ?nfi/(n  +  p)a]  (4.13) 

K2[lll]  =  ^-[(ni  +  3p)/4a  +  2{n  -  p)/4a  +  4nR/(«  +  p)a] 


A  knowledge  of  the  parameter  £  =  (a  -  B)/a  +  ft  is  enough  for  our 
calculations  and  so  individual  a  and  ft  are  not  evaluated.  As  has  been 
shown  in  Eq.  (4.4) 


The  parameters  can  be  expressed  in  terms  of  y ,  <?  and  c  as  follows: 

K3C100]  =  y  -  f  +  9r  (4.15) 

K3n  10]  -  3C12  f}  =  Dy  -  ES  +  Fe  (4.16) 

where 

D=  1  +3(7  -  r)2 
E  =  1  +  3(1  +  r)2 

and 

F  =  3[4  +  2 r,  -  f.2] 

K3 [ m ]  -  4C12  f}  =  Hy  -  P6  +  Qe  (4.17) 

where 

H  =  1  +  4(1  -  rj2  +  (1  -  O3 

p  =  l  +  j  (1  +  o2 

0=l+^-(f2  +  6f;  +  5)  . 


Solving  between  (4.15),  (4.16)  and  (4.17)  we  get  expressions  for  r,  f  and 
v  as  follows: 
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where 


G  =  K3[H0]  -  3C-| 2 

and 

R  =  K3[111]  -  4C-j  2  62 


(k3[ioo]  •  d  -  g)  (9D  _  F) 

6  =  Cr_:~D]  (E  -  0)  c 

and 


y  =  K3[100]  +  6  -  9e  . 


(4. IP) 


(4.19) 


(4.20) 


Substituting  these  force  constants  in  Eq.  (4.10)  leads  to  the  six 
independent  TOE  constants  for  diamond-like  solids.  Evaluation  of  r, 
y  and  6  as  a  function  of  temperature  results  in  the  TOE  constants  as  a 
function  of  temperature. 

c.  Results  -  Temperature  Variation  of  the  TOE  Constants  of  Silicon 
and  Germanium 

As  described  in  Chapter  III,  we  have  measured  the  nonlinearity 

parameters  and  hence  the  «3  parameters  along  the  three  principal  directions 

of  silicon  as  a  function  of  temperature  from  room  temperature  to  3°K 

using  the  ultrasonic  harmonic  neneration  technioue.  The  resul ts  are  uiven  in 

Chapter  III.  Using  the  same  technique, the  nonlinearity  parameters  and  K3 

46  59 

parameters  of  germanium  have  been  measured  by  Rains  ’  in  the  same 
temperature  range.  The  SOE  constants  of  silicon  and  germanium  have  also 


t 


d 
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been  measured  in  the  entire  temperature  region.  From  these  measured 
and  K3  values,  substituting  in  (4.14),  (4.18),  (a.19)  and  (4.20)  we 
have  evaluated  the  force  constants  £,  e,  y  and  s  as  a  function  of 
temperature  for  both  silicon  and  germanium.  In  Table  8  we  have  tabulated 
the  results. 

In  Figures  23  and  24  we  have  drawn  the  temperature  variation  of 
the  anharmonic  force  constants  of  silicon  and  oermanium.  It  may  be 
noted  from  the  figures  that  the  force  constants  are  susceptible  to  con¬ 
siderable  temperature  variations  at  low  temperatures  assuminn  nearly 
constant  values  at  higher  temperatures.  The  harmonic  force  constants  do 
not  vary  much  with  temperature.  The  quantity  r  which  is  defined  bv 
(4.7)  appearing  in  the  set  of  equations  (4.10)  is  almost  indeDendent  of 
temperature,  as  can  be  seen  from  Table  8. 

The  force  constants  at  various  temperatures  have  been  substituted 
in  Eqs.  (4.10)  to  evaluate  all  the  six  independent  TOE  constants  of 
silicon  and  germanium  as  a  function  of  temperature  between  300  and  3°K. 

In  Table  9  we  tabulate  the  results.  In  Figure  25  we  draw  the  TOE 


constants  and  C-jgg  of  silicon  as  a  function  of  temnerature. 

Temperature  variation  of  -^3 »  ^44  and  C455  of  silicon  are  drawn  in 

Figure  26.  Similarly,  in  Figures  27  and  28  we  have  drawn  the  TOE 
constants  of  germanium  as  a  function  of  temperature. 

In  the  figures  we  have  plotted  both  calculated  data  points  and 
a  best  fit  curve  through  them  to  show  the  effect  of  error  propagation. 


Even  in  the  curves  showing  greatest  scatter  of  data,  namely  the 
curves,  the  temperature  dependence  seems  very  clear.  As  can  be  seen 
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Table  8.  Temperature  Variation  of  the  Force  Constants  of  Silicon  and 
Germanium 


Temperature 

°K 

2C12 

x  1012 

x  io'2 

-Y 

ln12 
x  1 0 

C11  +  C12 

Silicon 

4.02 

0.5594 

0.5630 

0.6273 

3.1061 

10.57 

0.5594 

0.5327 

0.5614 

3.4246 

16.14 

0.5593 

0.5174 

0.5510 

3.5320 

24.93 

0.5593 

0.5089 

0.5246 

3.6228 

35.08 

0.5593 

0.4940 

0.4856 

3.7813 

43.49 

0.5592 

0.4926 

0.4557 

3.8137 

51.78 

0.5592 

0.4742 

0.4420 

4.0031 

61.04 

0.5592 

0.4434 

0.4192 

4.1697 

77.61 

0.5591 

0.4342 

0.^220 

4.1629 

100.42 

0.5589 

0.4364 

0.3998 

4.1511 

118.77 

0.5588 

0.4378 

0.3927 

4 . 1 344 

131.51 

0.5586 

0.4343 

0.3877 

4.1581 

147.89 

0.5585 

0.4396 

0.3840 

4.1075 

167.99 

0.5582 

0.4301 

0.3760 

4.1874 

197.93 

0.5579 

0.4254 

0.3608 

4.2228 

223.37 

0.5575 

0.4268 

0.3536 

4.1992 

241.70 

0.5573 

0.4233 

0.3474 

4.2214 

260.13 

0.5571 

0.4247 

0.3452 

4.1969 

272.48 

0.5570 

0.4225 

0.3392 

4.2137 

Room  temp. 

0.5567 

0.4244 

0.3357 

4.1800 

Germanium 


3.0 

0.5476 

0.6893 

0.61 5^ 

1 .6509 

9.0 

0.5475 

0.5985 

0.4839 

2.5996 

19.3 

0.5472 

0.5401 

0.4169 

3.1922 

27.8 

0.5471 

0.5494 

0.3637 

3,1617 

43.2 

0.5469 

0.5476 

0.2699 

3,2617 
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0.5464 

0.5932 
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2,9026 
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0.4776 
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100.5 
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0.5459 
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0.5458 

0.5319 

0.1948 

3,3781 

184.5 

0.5458 

0.4550 

0.1622 

4.0228 

196.5 

0.5458 

0.5029 

0.1948 

3.5791 

220.5 

0.5457 

0.4968 

0,1896 

3.6392 

245.1 

0.5455 

0.4522 

0.2276 

3,8726 

270.2 

0.5454 

0,3221 

0.2875 

4.6936 

Room  temp. 

0.5459 

0.3496 
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Figure  24.  Temperature  variation  of  the  Keating  anharmonic  force 
constants  of  germanium. 
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Figure  25.  Temperature  variation  of  the  TOE  constants 
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from  the  figures,  the  temperature  variation  of  the  TOE  constants  of 
silicon  and  germanium  are  similar,  as  can  be  expected.  The  room 
temperature  values  are  in  agreement  with  those  obtained  bv  McSkimin 

pc 

et  al.  The  greatest  disagreement  is  found  in  C^.  The  constants 
which  show  maximum  sensitivity  to  temperature  chanoes  are  C-^  and  C-^. 
They  take  on  positive  values  at  very  low  temperatures  for  both 
materials . 

It  should  be  mentioned  that  in  the  above  calculations  one 
obviously  has  to  make  the  assumption  that  the  thermal  effects  in  elastic 
constants  are  solely  due  to  changes  in  force  constants.  Th i c  assumption 
is  open  to  question.  Temperature  dependent  terms  in  the  internal 
energy  also  might  contribute  to  the  elastic  constants.  Such  effects,  we 
believe,  would  be  negligible  compared  to  the  variations  occurring  from 
changes  in  force  constants. 

d_.  Temperature  Dependence  of  the  Griineisen  Parameter,  Thermal 
E x pa n s_i on  and  Other  Anharmonic  Parameters  of  Si  1 i con  a nd 
Germanium 

As  is  well  known,  the  Griineisen  parameter  represents  the  strain 
derivative  of  the  lattice  vibrational  frequencies  and  is  an  important 
quantity  in  describing  the  anharmonic  properties  of  solids.  The  method 
of  evaluating  the  Griineisen  parameter  usinq  quasiharmonic  theory  is 

71  72 

well  established  in  literature.  ’  The  Gruneisen  parameter  is 
defined  by 


KfCv  Kc-Cp 


(4 . ?0) 


where  <  is  the  thermal  expansivity,  Ky  and  K<-  are  the  isothermal  and 
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adiabatic  comnressi bi 1 i ties  and  Cp  and  Cy  are  the  isochoric  and  isobaric 

heat  capacities.  Generalized  y ' s  measure  the  strain  dependence  of  the 

phonon  frequencies  and  are  no  lonoer  freouenev  dependent  in  the  Debve 

model  where  the  lattice  vibrations  are  replaced  by  standina  wave  modes 

of  a  dispersionless  elastic  continuum.  Their  relationship  with  SOF  and 

58 

TOE  coefficients  have  been  derived  earlier  by  Bruqger  for  arbitrary 
crystal  symmetry  and  have  been  specialized  to  cubic  point  groups. 

In  the  quasiharmonic  approximation,  the  Gruneisen  y  can  also  be  expressed 
as  the  weighted  average  of  the  generalized  Gruneisen  parameters  ■ ■  by  the 
relation 

7  =  a  YiCi/)Ci  (4.21) 

where  expresses  the  volume  dependence  of  the  lattice  vibrational 
frequency  for  the  ith  mode  and  C.  is  the  Einstein  heat  caoacity  associated 
with  that  mode.  If  the  sound  speeds  and  their  stress  derivatives  or  the 
SOE  and  TOE  constants  are  known,  Eq.  (4.21)  can  be  evaluated  in  the 
continuum  model  even  if  the  dispersion  curves  alona  many  directions  in 
the  stressed  and  unstressed  crystal  are  not  known. 

72 

We  follow  the  method  due  to  Brugper  and  Fritz  to  determine  the 
Gruneisen y  '  s  and  thermal  expansion  of  silicon  and  germanium.  If  p  stands 
for  the  branch  index  and  o  *or  the  wave  vector,  Eq.  (4.21)  can  be 
written  in  the  form 

=  /  y(p»q)C(n,a)/  7  C(p,a) 

Pd  pq 


Y 


(4,2?) 
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The  continuum  model  has  been  introduced  for  the  evaluation  of  v  with 
the  following  assumptions: 

(i)  The  excitation  of  the  oDtic  modes  may  be  neglected,  so  that 
the  branch  index  p  takes  only  the  values  1,  2  and  3. 

(ii)  The  acoustic  modes  obey  the  Debye  distribution  function 
(per  unit  volume)  qiven  by 

g(p,q)dqd  =  (— )q2dqdi?  .  .22) 

(iii)  The  maximum  value  of  the  wave  vector  o  along  any  direction 
is  equal  to  the  Debye  radius 

qD  =  (6v~)1/3  (4.24) 

where  is  the  volume  of  the  primitive  cell. 

(iv)  The  acoustic  modes  are  either  nondispersive  or  they  obey 
the  sinusoidal  dispersion  relation 


ir  Q 


•  t » 

=  sin(-?-  - 


(A. 25) 


max  ~  “'D 

of  the  Born-von  Karman  model . 

(v)  The  generalized  Gruneisen  parameter  is  independent  of  the 
wave  number  and  is  aiven  bv  their  long  wave  limits. 

Applying  the  above  assumptions,  Eq.  (4.22)  takes  the  form 


Y  =  n  dMP,N)C(P,N)/E  OdnC(P.N) 


(4.26) 


wi  th 


qn  3  ^ 

C(P,N)  =  (^)J 


dr. 


0 


rYeQ 

(e°-l)? 


(4.27) 
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is  the  Debye  radius  and  in  the  Debye  model  the  function  0  is  given 
by 

0(P,N,O  =  [o(P,Nl/T]f.  (4.28) 

where  f.(P,N)  is  the  Debye  characteristic  temperature  of  the  pth  mode 

along  the  direction  of  q  specified  by  the  unit  vector  N.  If  S(P,N)  is 

78 

the  elastic  wave  velocity  of  the  pth  mode,  then 

e(P.N)  =  (^)S(P,N)  .  (4.29) 

For  cubic  crystals  with  a  fourfold  axis  of  symmetry,  the 

72 

expressions  for  y(P,N)  have  been  derived  earlier  and  are  given  by 

y(P,N)  =  -  ~  [3B  +  2W  +  K]  (4.30) 

where 

W(P,N)  =  C]1K1  +  C44K2  +  C12K3  (4.31) 

and 

K(P,N)  =  CjKj  +  C2K2  +  C3K3  .  (4.32) 

The  unknown  quantities  appearing  in  (4.31)  and  (4.32)  are 

«i  (P,N)  =  N^Uj2  +  U22U22  +  N32U32  | 

K2(P,N)  =  (N2U3  +  N3U2)2  +  (N3U1  +  N1U3)2  +  (N,U2  +  N2U1)2|(4.33) 

k3(p,n)  =  2(n2n3u?u3  +  n1n3u]u3  +  n1n2u1u?  ] 


and 


ci  =  cm  +  2cii2 
C2  =  C 1 44  +  2C166 
C3  =  C1 23  +  2C1 1 2 


(4.34) 


C..'s  and  C...  ,  are  the  SOE  and  TOE  constants  of  the  solid.  Also 

1  J  T  J  K  S 


B  r  3^cn  +  2ci 2 ' 


(4.35) 


is  the  bulk  modulus.  The  N's  and  U's  appearing  in  Eqs.  (4.33)  are  the 
direction  cosines  for  the  direction  of  propagation  and  direction  of 
polarization  characterized  by  the  branch  index  P. 

The  thermal  expansion  is  given  by 


^  =  K  l  y(p,q)C(p,q) 

pq 


(4.36) 


where  K  is  the  isothermal  compressibility.  The  thermal  exoansion  of 
cubic  crystals  is  isotropic  and  hence  only  the  scalar  y  of  Eq.  (4.26) 
need  be  evaluated  for  its  determination. 

Followinq  the  above  procedure  we  have  determined  the  Gruneisen 
parameter  of  silicon  and  germanium  as  a  function  of  temperature  from 
300  to  3'K  usinq  the  temperature  dependent  SOE  and  TOE  constants  that 
we  have  determined  in  the  previous  section.  In  Fioures  29  and  30  we 
plot  the  Gruneisen  parameter  results  obtained  as  a  function  of  temperature 
for  silicon  and  germanium  respectively.  The  computations  have  been 
done  with  the  help  of  a  computer.  Our  curves  are  designated  as  "present 
work"  in  the  figures.  With  the  same  computer  program  we  have  also 
evaluated  the  temperature  variation  of  y  usinq  room  temperature  values 
of  TOE  constants.  This  is  the  method  adopted  by  Bruqger  and  Fritz7'  who 
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Figure  29.  Griineisen  parameter  of  silicon  plotted  as  a  function 
of  temperature. 
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have  assumed  that  TOE  constants  are  independent  of  temperature.  These 

curves  designated  as  "Bruqoer  and  Fritz  Method"  are  also  given  in  the 

figures  along  with  the  original  curves  drawn  by  Pruager  and  Fritz 

themselves  for  comparison.  The  curves  drawn  by  Brugqer  and  Fritz  are 

denoted  by  "Brugger  and  Fritz  Curves."  The  small  difference  between  our 

curves  drawn  by  the  Brugger  and  Fritz  method  and  the  curves  drawn  by 

Brugger  and  Fritz  are  due  to  slight  difference  in  the  numerical  values 

used.  We  have  used  the  room  temperature  values  of  the  TOE  constants 

26 

reported  by  McSkimin  et  al .  for  both  silicon  and  qermanium,  whereas 

Brugger  and  Fritz  used  the  pressure  derivatives  of  the  SOE  constants 

54  79 

reported  by  McSkimin  et  al.'  for  silicon  and  McSkimin  et  al.  for 

germanium.  In  the  figures  we  have  also  drawn  the  experimental  curves 

69 

of  Y  obtained  from  thermal  expansion  measurements.  These  curves  are 
denoted  by  "Thermal  Exp"  in  the  curves.  It  can  be  seen  from  the  curves 
that  our  curves  which  take  into  account  the  temperature  variation  of 
the  TOE  constants  are  in  better  agreement  with  experimental  curves, 
even  though  they  also  do  not  show  any  negat^'e  values  for  the  Griineisen 
parameter  which  has  been  calculated  from  thermal  expansion  data.  A 
realistic  physical  reasoning  for  this  disagreement  is  included  in  the 
discussion  section. 

In  the  quasi  harmonic  approximation,  the  thermal  expansion  is 
proportional  to  the  Griineisen  parameter  (Eg.  4.36)  and  so  thermal 
expansion  follows  a  similar  temperature  variation  as  the  Griineisen 
parameter.  Since  the  thermal  expansion  curves  are  similar  to  the 
Griineisen  parameter  curves,  they  are  not  reproduced.  In  the 
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high  temperature  limit  the  value  of  y  calculated  from  thermal  expansion 
agrees  well  with  measured  values. 

We  have  determined  and  drawn  as  a  function  of  temperature  the 
isothermal  first  pressure  derivatives  of  the  adiabatic  SOE  constants 
from  the  temperature  dependent  TOE  constants.  The  relationships 
between  pressure  derivatives  of  the  SOE  constants  and  TOE  constants  are 

i i  i  80 ,81 

well  known 


dC 


11 


dP 


2Cn  +  2Ci2  +  cm  +  2cn  2 , 
(  C11  +  2C1 2  "  * 


dCl 2  _  ,'C11  '  C1 2  +  C 1 23  +  2C1 1 2 x 

l  r  a.  or  ' 


dP 


dC 


44 
dP  ‘ 


C11  +  2C12 


, Cn  *  2C1 2  +  C44  +  C 1 44  +  2C166, 
CH  +  2C12 


(4.37) 


dB  _  /Cm  +  6C112  +  2C1 23  N 
dP  “  '  1  3(C11  +  2C12)  ; 

where  B  =  2(C-|  -|  +  2C ^ 2 )  is  the  bulk  modulus.  The  above  pressure 

derivatives  drawn  as  a  function  of  temperature  for  silicon  and  germanium 

are  reproduced  in  figures  31  and  32  respectively. 
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Anderson  stated  that  the  adiabatic  bulk  modulus  of  a  material 
varies  with  temperature  according  to 


BS  B0  "  bT  e 


VT 


(4.38) 


where  B<.  is  the  bulk  modulus  at  temperature  T  and  B^  is  the  bulk  modulus 
at  absolute  zero,  b  and  T^  are  empirical  parameters  characteristic  of 
the  solid.  Anderson  showed  that  the  constant  b  is  given  by 
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where  is  the  specific  volume  per  atom  at  absolute  zero,  I'  is  the  gas 

constant,  r  is  the  Gruneisen  parameter  and  1  is  a  constant  .ailed  the 

Anderson-Gruneisen  parameter  and  it  plays  an  imnortant  role  in  descrihine 

the  temperature  denendence  of  the  bulk  modulus  and  related  anharmonic 
83 

properties.  Chanq  showed  that  4  is  related  to  the  pressure  derivative 
of  the  bulk  modulus  as 


hp  *  ^ 


[4.40) 


Usinq  this  definition,  the  relationship  between  ■  and  TOE  constants  has 
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been  obtained  by  Rao  as 


This  expression  has  been  used  to  determine  the  temperature  variation  of 
■  usinq  temperature  dependent  TOE  constants.  r  as  a  function  of  tem¬ 
perature  for  silicon  and  germanium  are  also  drawn  in  Figures  31  and  3?. 


e.  Discussion  and_ Conclusion 

An  elegant  method  of  combining  an  experimental  result  and  a 
successful  theoretical  model  has  enabled  us  to  determine  the  temperature 
dependence  of  the  TOE  constants  of  silicon  and  germanium  from  liquid 
helium  to  room  temperature  for  the  first  time.  It  is  a  remarkable 
step  towards  understanding  the  low  temperature  anharmonic  properties 
of  these  solids  even  thouoh  the  results  do  not  fully  account  for  some 
of  the  anomalous  thermal  properties  exhibited  by  them  at  low  temperatures. 
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Our  assumption  that  the  temperature  variation  of  the  elastic  constants 
are  due  to  thermal  effects  in  the  force  constants  is  justifiable  as  we 
see  that  our  f.rune i sen  parameter  curves  aoree  with  the  curves  of  others. 

As  a  matter  of  fact,  our  curves  are  in  better  aoreement  with  curves 
derived  from  experimental  thermal  expansion  curves  than  any  nrevious 
resul ts . 

The  Gruneisen  parameter  and  thermal  expansion  of  silicon  and 
qermanium  are  nenative  at  certain  temperatures  in  the  low  temperature 
reqion.  This  anomalous  property  has  attracted  a  nreat  amount  of  attention. 
But  none  of  the  calculations  based  on  the  anisotropic  continuum  model 
could  predict  or  fully  account  for  this  anomalous  behavior.  Now  it 
is  more  or  less  well  established  that  the  negative  thermal  properties 
are  a  result  of  some  of  ?he  mode  r»runef$pn  parameters  ooinq  nenative  or 
the  phonon  frequencies  decreasinn  with  stress.  A  straightforward 
calculation  based  on  anisotropic  continuum  model  in  the  quasi  harmonic 
approximation  cannot  account  for  these  negative  mode  ' s .  The  ove'all 
temperature  dependence  o*  ■  is  explained  as  follows:  At  very  low 
tempera  turns ,  the  region  in  which  .  is  positive,  the  elastic  portion 
of  the  lattice  spectrum  dominates  uivinn  positive  ..  As  the  temperature 
is  increased,  the  entire  transverse  acoustic  modes  (TA)  contribute 
appreciably  to  v.  The  mode  . 's  for  the  transverse  modes  decrease  with 
temperature.  In  the  negative  •  temperature  reqion,  the  contribution 
from  TA  modes  dominates  over  the  ocher  modes  making  the  net  effect, 
negative.  As  the  temperature  is  increased  further,  the  remainder  of  the 
acoustic  branches  coupled  with  the  optic  modes  make  .  positive  again.  In 


HO 


the  high  temperature  region  all  the  theoretical  calculations  agree 

with  experimental  results.  A  complete  understanding  of  the  situation 

demands  complete  experimental  data  on  the  volume  dependence  of  the 

frequencies  of  the  individual  branches  in  the  dispersive  region  in  the 

entire  temperature  region.  Since  such  data  are  not  available, 
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Bienenstock  has  done  a  calculation  for  germanium  based  on  a  modified 
she  i  model  using  the  volume  dependence  of  the  elastic  moduli  and 

OC 

phonon  frequencies  reported  by  Brockhouse  and  Iyengar.  Bienenstock 

introduces  an  induced  dipole  term  which  varies  as  V~n  where  n  is  an 

adjustable  parameter  nearly  egual  to  1  and  includes  the  optic  modes  in 

the  Einstein  approximation.  He  has  calculated  the  Gruneisen  parameters 

y  ■  for  a  number  of  acoustic  modes  and  has  found  that  the  transverse 

modes  are  rapidly  varying  functions  of  the  wave  vector  q.  They  are 

assumed  to  be  positive  for  small  o  and  negative  near  the  Brillouin  zone 

surface.  The  calculated  y  is  in  agreement  with  experiment  at  all 

temperatures  when  ontical  modes  are  included.  Studies  on  phonon 

assisted  tunnelling  under  pressure  and  measurements  of  the  stress- 

induced  shifts  of  the  threshold  energy  for  tunnellinn  in  germanium  and 

subsequent  calculation  of  for  TA  modes  near  the  the  Brillouin  zone 
87 

boundary  by  Pyne  confirms  the  existence  of  negative  y.  for  TA  modes. 

General  Conclusion 

The  work  presented  in  this  technical  report  forms  a  nearly 
complete  study  of  the  temperature  dependence  of  the  third-order  elastic 
constants  of  the  typical  diamond-like  solids  silicon  and  germanium. 


The  harmonic  generation  technique  has  found  new  application  to  the 
study  of  the  anharmonic  properties  of  these  solids.  Similar  works 
applied  to  other  groups  of  solids  might  lead  to  an  understanding  of 
the  temperature  variation  of  their  third-order  elastic  constants  and 
other  anharmonic  properties. 
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APPENDIX  A . 1 


THE  DETECTOR  EQUIVALENT  CIRCUIT 

Figure  33  shows  the  detector  as  it  is  connected  for  measurements 
of  nonlinearity.  It  is  assumed  that  the  detector  load  is  an  ideal 
resistive  load  R.  A  larqe  coupling  capacitor  C^  provides  dc  blocking 
for  the  amplifier.  Its  effect  on  the  RF  sionals  can  be  ignored.  is 
a  large  (%  1  to  10  mb)  resistor  used  to  prevent  larqe  current  discharges 
in  the  event  of  arcing  in  the  detector.  It  is  very  large  compared  to 
R  so  that  no  appreciable  current  flows  through  it  and  V^.  If  one 
considers  the  sample  end  vibrating  sinusoidally  at  frequency  c>  with 
amplitude  2A,  where  A  is  the  acoustic  wave  amplitude  in  the  sample, 
then  the  gap  spacing  chanqes  with  time  according  to 

S  =  SQ  +  2A  sin  mt  .  (A. 1-1) 

The  factor  2  enters  because  the  reflection  of  the  wave  at  the  stress- 
free  surface  leads  to  a  vibration  amplitude  of  the  sample  end  which  is 
twice  the  wave  amplitude  inside  the  sample.  Thus  the  capacitance  of  the 
detector,  considered  as  a  parallel  plate  capacitor,  is  oiven  by 

1/C  =  1/Cn(l  +  |A  sin  ...t)  (A.  1  -?) 

u 

where  Cn  =  ^  is  the  static  capacitance  of  the  detector,  .  ,  .  and  Sn 
u 

being  the  dielectric  constant  of  the  medium  (=1  for  air),  area  of  the 
electrode  and  static  spacing  of  the  detector,  respectively.  Summing 
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voltaqes  around  the  outside  loop  of  Figure  33,  we  qet 

^  -  vb  =  -R  af  =  v(t)  .  (A. 1-3) 

Substituting  ( A .  1  - 1 ) and  (A. 1-2)  in  (A. 1-3)  and  assuming  a  steady-state 
solution  of  the  form 


(A. 1-4) 


(A. 1-5) 


f  =  tan'1  (R  ...  CQ)  .  (A. 1-6) 

2AV 

Consider  a  signal  generator  of  amplitude  V  =  feedinq  R 

0 

through  a  capacitance  of  value  CQ.  The  voltage  across  R  can  be  shown 

to  be  precisely  the  same  (in  both  amplitude  and  phase)  as  Eg.  (A. 1-5).  A 

complex  load  of  impedance  Z,  rather  than  R,  can  also  be  treated.  One  has 
AVb 

simply  sim.it  feedinq  Cn  and  Z  in  series.  This  situation  more 
^0  u 

accurately  represents  the  experimental  case  because  of  stray  line 
capacitance  shunting  the  input  to  the  amplifier.  For  most  experimental 
situations  this  stray  capacitance  is  sufficiently  large  to  reouire  that 
during  the  calibration  to  determine  V,  one  feeds  through  a  substitutional 
capacitor  of  value  C^,  even  if  R  is  chosen  to  be  quite  large. 


APPENDIX  A. 2 


RF  GATE,  MOSFET  GATE  AND  FREOUENCY  DOUBLER 

To  prevent  overloading  of  the  amplifiers  due  to  the  feedthrouqh 
signal,  two  gates  have  been  built.  An  RF  gate,  shown  in  Fiaure  3fl,  is 
used  between  the  capacitive  detector  and  the  30  MHz  amplifier  to  pass 
only  the  echoes  of  interest.  A  slower  MOSFET  gate,  shown  in  Fiqure  35, 
is  used  between  the  detected  output  of  the  30  MHz  or  60  MHz  broadband 
amplifier  and  the  boxcar  integrator.  The  RF  gate  is  able  to  handle  the 
RF  pulses,  but  the  on-to-off  ratio  is  not  very  great.  However,  if  a 
xlO  attenuator  is  placed  between  the  gate  and  the  30  MHz  amplifier, 
and  the  amplifier  is  operated  at  higher  gain  (to  prevent  overloading  of 
the  first  stage,  which  has  fixed  gain),  the  attentuation  of  the  feed¬ 
through  is  sufficient  to  prevent  overload  of  the  30  MHz  amplifier. 

The  MOSFET  gate  is  used  to  protect  the  more  sensitive  boxcar 
integrator  from  overload.  The  on-to-off  ratio  of  the  MOSFET  gate  is 
very  hiqh,  but  the  gate  is  useful  only  for  low  frequency  signals.  Both 
gates  are  driven  bv  a  pair  of  complementary  pulses  from  a  GR1398-A  pulse 
generator.  A  delayed  triqqer  output  from  the  boxcar  integrator  allows 
the  gatinq  pulse  to  be  positioned  in  time. 

The  frequency  doubler  used  in  the  calibration  equipment  during 
room  temperature  measurements  is  essentially  a  ring  bridge  modulator- 
phase  sensitive  detector  followed  by  a  tunable  series  resonant  filter 
shown  in  Fiqure  36.  It  provides  a  signal  which  is  precisely  twice  the 
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frequency  of  the  fundamental  used  in  the  study.  The  diodes  D ^ - D ^  shown 
in  the  fiqure  are  RCA  type  CA3039.  The  doubler  has  been  tested  for 
spectral  purity  with  a  spectrum  analyzer. 


APPENDIX  A.  3 


ATTENUATION  AND  HARMONIC  LOSS  EFFECTS 

As  the  ultrasonic  wave  propagates  in  the  sample,  the  second 
harmonic  amplitude  A2  chanqes  as  a  function  of  distance.  There  is  a 
growth  of  A2  proportional  to  the  square  of  the  first  harmonic  amplitude 
A.|  and  there  is  a  decay  of  the  wave  due  to  attenuation  of  Aj.  Therefore 
the  amplitude  of  the  second  harmonic  as  a  function  of  distance  in  the 
sample  a  is  given  by  the  differential  eouation 


dA2 

da 


*2A2  +  §  k 


(A. 3-1) 


where  ^  i 5  the  attenuation  coefficient  of  A,,.  If  one  assumes  that  the 
change  in  A^  as  a  function  of  a  is  due  only  to  the  attenuation  of  A-j  and 
the  loss  of  power  into  the  second  harmonic,  neglecting  the  power  loss 
into  third  and  higher  harmonics,  one  obtains 


dA-j 

da 


-<1A1 


2  k2 

8  k 


A?A1 


(A. 3-2) 


where  u-j  is  the  attenuation  coefficient  of  A-j  . 

The  simultaneous  solution  of  Eqs.  (A. 3-1)  and  (A. 3-2)  has  not 
been  obtained,  but  two  special  cases  can  be  solved.  The  first  special 
case  is  the  one  for  which  we  can  neglect  the  term  due  to  power  loss  from 
the  fundamental  into  the  second  harmonic.  Equation  (A. 3-2)  then  becomes 


dA-j 

Ha 


'rtlAl 


(A. 3-3) 


156 


157 


Equations  (A. 3-1)  and  (A. 3-3)  can  be  solved  simultaneously  to  give 


A 


2 


2' 1 1  a 
e 


-  2  1 1  a  -a 

e  -  e 


*2 


(A. 3-4) 
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The  attenuation  coefficients  for  silicon  have  been  measured. 

For  a  30  MHz  fundamental,  the  coefficients  are  of  the  order  of 

i ,  £  0.09  db/cm  and  £  0.27  db/cm.  For  a  5  cm  sample  (which  is  lonqer 

2 

than  any  of  the  samples  used  in  the  present  measurements)  A^/A^  qives 
2 

a  value  for  ( 3/8 ) k  Fa  which  is  approximately  2.6"  low.  Since  this  is 
much  smaller  than  the  other  sources  of  error  in  the  absolute  measurements, 
no  attenuation  correction  is  necessary. 

The  second  special  case  is  that  for  which  =  0. 

Equations  (A. 3-1)  and  (A. 3-2)  then  reduce  to 


^2  32.2 

~da”  "  8  k  '  A1 


(A. 3-5) 


and 


dA  i  n  2 

da  8  1  M2  1 


(A. 3-6) 


Equations  (A. 3-5)  and  (A. 3-6)  have  the  simultaneous  solution 


2  2  sin  h'^ (A^/A^ ) 

£-  k  i-:a  =  - — — - — 

8  /A— 2 - 2 

V  h i  +  a2 


i-V'a  ^  +  "  2' 


en  (A2  +^A2  +  A1  )  -  v.n  A 


(A. 3-7) 
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This  equation  is  the  same  as  Eq.  (3.2)  of  Chapter  III.  For  the  entire 

range  of  amplitudes  used  in  the  room  temperature  measurements  of 

2 

silicon,  the  difference  between  ( A^/A^  )  and  the  quantity  calculated 
from  Eq.  (A. 3-7)  is  very  small. 
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